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Abstract 

Artificial interface conditions parametrized by a complex number 80 are introduced for 
1D-Schr6dinger operators. When this complex parameter equals the parameter 9 £ iR of the 
complex deformation which unveils the shape resonances, the Hamiltonian becomes dissipative. 
This makes possible an adiabatic theory for the time evolution of resonant states for arbitrarily 
large time scales. The effect of the artificial interface conditions on the important stationary 
quantities involved in quantum transport models is also checked to be as small as wanted, in 
the polynomial scale (h )jyeM as h — > 0, according to 8q. 



1 Introduction 

The adiabatic evolution of resonances is an old problem which has received various answers in the 
last twenty years. The most effective results where obtained by remaining on the real spectrum 
and by considering the evolution of quasi- resonant states (see for example [4"8][60"][T]). Motivated 
by nonlinear problems coming from the modelling of quantum electronic transport, we reconsider 
this problem and propose a new approach which rely on a modification of the initial kinetic en- 
ergy operator — h 2 A into —h 2 Ag g where 0q parametrizes artificial interface conditions. With this 
analysis, we aim at developing reduced models for the nonlinear dynamics of transverse quantum 
transport in resonant tunneling diodes or possibly more complex structures. A functional frame- 
work for such a model has been proposed in [45] and implements a dynamically nonlinear version of 
the Landauer-Biittikcr approach based on Mourre's theory and Sigal-Soffer propagation estimates 
(see [TB] and [TS] for an alternative presentation of the stationary problem). The derivation of 
reduced models for the steady state problem has been developed in [20 [2TJ [IS] on the basis of 
Hclffer-Sjostrand analysis of resonances in [30]. This asymptotic model elucidated the influence of 
the geometry of the potential on the feasibility of hysteresis phenomena already studied in [33] [50] . 
Numerical applications have been carried out in realistic Ga — As or Si — S1O2 structure in |19) 
and [18] . showing a good agreement with previous numerical simulations in |17) or [41] and finally 
predicting the possibility of exotic bifurcation diagrams. From the modelling point of view a diffi- 
culty comes from the phase-space description of the tunnel effect which can be summarized with 
the importance in the asymptotic nonlinear system of the asymptotic value of the branching ratio 

K^Vi(+fc,.),<^)| 2 

U = nm ; . (1.1) 

In the above formula zj 1 = — iT^ is a resonance for the Hamiltonian H h = —h 2 A + V — W h with 
a semiclassical island V and a quantum well W h , ip-(±k, .) are the generalized eigenfunctions for 
the filled well Hamiltonian H" = —h 2 A + V with a momentum ±fc such that k 2 ~ Ej and $^ is 
the j-th eigenfunction of the Dirichlet Hamiltonian = —h 2 A + V — W h on some finite interval 
(a, b). The imaginary part of the resonance is given by the Fermi Golden rule proved in |21j 

I (V^V 1 (+£;,•),</>?) I 2 \(W h ^ h (-k,-),^)\ 2 

r , (1 + o(1)) = f^iM + , (1.2) 
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where the interaction with the continuous spectrum leading to a resonance contains two contribu- 
tions from the left-hand side with +k and from the right-hand side with — k. 

Our purpose is the derivation of reduced models for the dynamics of quantum nonlinear systems 
like it has been done in the stationary case in [20] [S] [M] with the following motto: The (nonlinear) 
phenomena are governed by a finite number of resonant states. As this was already explained in 
[20] , such a remark dictates the scaling of the potential which leads to the small parameter analysis 
and in the end to effective reduced models even when h ~ 0.1 or 0.3 (see [T9]). 
For the dynamical problem, the time evolution of resonant states have to be considered possibly 
with a time-dependent potential. And it is known that this is a rather subtle point. Quantum 
resonances follow almost but not exactly the general intuition (see [53) 1 and remain an inexhaustible 
playground for mathematical analysis. For example, the exponential decay law is an approximation 
which has a physical interpretation in terms of the evolution of quasi-resonant states (truncated 
resonant states which lie on the real L 2 -space) and writes as 

where the remainder term R(t) is small only in the range of times scaled as jr-. A very accurate 
analysis of this has been done in 27 [57 58 [59 39 [38 and adiabatic results for slowly varying 
potentials and for quasi-resonant states have been obtained in [38] (see also [50] and pQ in a similar 
spirit) on this range of time-scales. On the other side the relation 

e -«* fc (»ty 3 .= e -« B *e-* r ^ 

holds without remainder terms when 9 £ iM. + parametrizes a complex deformation of H according 
to the general approach to resonances (see [3][Tn][8T][23 [30 29 ). However, and this is well known 
within the analysis of resonances, the deformed generator iH h {9) is not maximal accretive although 
its spectrum lies in {Re z > 0} and no uniform estimates are available on e~ ttH ^ . One of the 
two next strategies have to be chosen: 

• Stay on the real space with quasi-resonant states, with uniform estimates of the semigroups, 
groups or dynamical systems (they are unitary) but with remainder term which can be 
neglected only on some parameter dependent range of time. 

• Consider the complex deformed situation and try to solve or bypass the defect of accretivity. 

Because the remainder terms seemed hard to handle within the original nonlinear problem and 
also because there may be multiple time scales to handle, due to the nonlinearity or due to several 
resonances involved in the nonlinear process, we chose the second one. 



In a one dimensional problem the simplest approach is the complex dilation method according 
to [5] |10] [23] . Since the possibly nonlinear potential with a compact support has a limited regularity 
inside the interval (a, b), this deformation is done only outside this interval following an approach 
already presented in [54] . The dilation is defined according to 

b) +b), x>b 

xe(a,b) (1.3) 
a) + a), x < a , 

and finally handled with 6 £ iK + . The conjugated Laplacian is 

with the domain made of functions u € H 2 (R \ {a, b}) with the interface conditions 

" e-iu(b+) = u{b-)- e-¥ u '(b+) = «'(&-) 

(1.4) 

e~zu(a~) = u(a + ); e~~u'(a~) — u / (a + ) . 



e^ip(e s (x - 
ip(x), 
e^ip(e e (x - 
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This can be viewed as a singular version of the black-box formalism of [55]. Additionally to the 
fact that this singular deformation is convenient for the original model with potential barriers 
presented with a discontinuous potential and with a nonlinear part inside (a, 6), the obstruction 
to the accretivity of Ug(— ft. 2 A)C/^ 1 is concentrated in two boundary terms at x = a and x = b in 



Re (it , iUg(- 



-h 2 A)U g - 1 u) m 



Re 



ih 2 (uu')\ _{e 



+ h 2 e 



-2 Re 



'sin(2Im 0) 



K\[o,6] 



\u'\ 2 dx. (1.5) 



Our strategy then relies on the introduction of artificial interface conditions parametrized with 9q 
which modify the operator — h 2 A. The parameter 0o is then chosen so that the above boundary 
term vanishes when 9 = 9q = it. The modified and deformed Hamiltonian Hg iT (9 = it) then 
generates a contraction semigroup and uniform estimates are available for e~ %tH i-r( %T ) Qr j or ^ ne 
dynamical system (U h (t,s)) t > s for time-dependent potentials. 

Hopefully this modification has a little effect on the Hamiltonian H h = — h 2 A + V — W h and all 
the quantities involved in the nonlinear problem, with explicit estimates with respect to 0o and h. 
Indeed all the quantities and even the exponentially small ones like Tj or the ones appearing in 
the branching ratio (jl.lj) experiment small relative variation with respect to 9q when 9o = ih N ° 
with N >5. 

In comparison with the modelling of artificial dissipative boundary conditions in [TT] [T3] [T5] [2] , our 
approach has the advantage of remaining close to the initial quantum model. Such a comparison 
is valuable and ensures the validity of numerical applications when the non-linear bifurcation 
phenomena are very sensitive to small variations of the data. 

Once the above comparison is done, it is checked that adiabatic evolution for a slowly varying 
potential or equivalently for the ^-dependent Cauchy problem 

ied t u = Hg o (0 O , t)u , u t=0 = u , 

with some exponentially large time scale | = e£, is adapted from the general approach for the 
adiabatic evolution of bound states of self-adjoint generators in [5] [15] [53]. 

Adiabatic dynamics have already been considered within the modelling of out-of-equilibrium quan- 
tum transport in [22 [7 9 playing with the continuous spectrum with self-adjoint techniques. Only 
partial results are known with non self-adjoint generators: in |42j only small time results are valid 
for resonances, in 44 bounded generators are considered and in [55] a general scheme for the 
the higher order construction of the adapted projector is done but without time propagation es- 
timates. In [55], A. Joye considered a general time-adiabatic evolution for semigroups in Banach 
spaces under a fixed gap condition and with analyticity assumptions: The exponential growth of 
the dynamical system ||S e (f, 0)|| < e~ is compensated by the 0{e~i) error of the adiabatic ap- 
proximation under analyticity assumptions (see 43 [34;) and lead to a total error 0( ) which 
is small when t < c/C. Our adaptation combines the uniform estimates due to the accretivity 
of the modified and deformed Hamiltonian with the accurate resolvent estimate provided by the 
accurate comparison with self-adjoint problems (shape resonances result from the coupling of some 
Dirichlet eigenvalues with a continuous spectrum). Although, the error associated with the adia- 
batic evolution is estimated at the first order as an <D(e 1 ~ 5 ), with S > as small as wanted, it is 
necessary to reconsider the higher-order method in |43] or |34j . because we work with small gaps 
(vanishing as h — > 0) and with non self-adjoint operators. Finally note that the exponential time 
scale is not necessarily related with the imaginary parts of resonances and several resonant states 
with various life-time scales are taken into account in our application. 



The outline of the article is the following. 

• The artificial interface conditions parametrized by 0o are introduced in Section [5] With 
these new interface conditions —A is transformed into a non self-adjoint operator conjugated 
with -A, M/(0 O )(-A)F4 7 (0 O )~ 1 with W{9 ) = ld L 2 +O(9 ). The case with a potential is 
illustrated with numerical computations. 
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• The functional analysis of the complex deformation parametrized with 9 is done in Section [31 
After introducing a Krein formula associated with the (9q, #)-dependent interface conditions, 
it mimics the standard approach to resonances summarized in [23] 29] but things have to be 
reconsidered for we start from an already non self-adjoint operator when 9$ 7^ 0. Assumptions 
on the time-dependent variations of the potential which ensure the well-posedness of the 
dynamical systems are specified in the end of this section. 

• The small parameter problem modelling quantum wells in a semiclassical island is introduced 
in Section |U Accurate exponential decay estimates are presented for the spectral problems 
reduced to (a, b) making use of the fact that our operators are proportional to — h 2 A outside 
[a, b]. 

• The Grushin problem leading to an accurate theory of resonances is presented in Section [5] 
There it is checked that the imaginary parts of the resonances, which are exponentially small, 
are little perturbed by the introduction of #o-dependent artificial conditions. The conclusion 
of this section is that all the quantities involved in the Fermi Golden Rule (|1.2p have little 
relative variations with respect to 9q, even the exponentially small ones. 

• Accurate parameter-dependent resolvent estimates for the whole space problem are done in 
Section [51 Again the Krein formula for the resolvent associated with the (#o, #)-dependent 
interface conditions is especially useful. 

• The adiabatic evolution of resonances is really introduced in [7J After specifying all the 
assumptions, the main result about this is stated in Theorem 17. II 

• The appendix contains various preliminary and sometimes well-known estimates, plus a vari- 
ation of the general adiabatic theory concerned with non self- adjoint maximal accretive op- 
erators in Section [Bl 



2 Artificial interface conditions 

Modified Hamiltonians with artificial interface conditions are introduced. It is checked that the 
effect of these artificial conditions parametrized by 6>o 6 C is of order \9q\ for all the quantities 
associated with the free Schrodingcr Hamiltonian —ft 2 A on the whole line. Instead of pursuing 
this analysis for general Schrodinger operators —h 2 A + V we simply give a numerical evidence of 
this stability with respect 9q . 



2.1 The modified Laplacian 

We consider a class of singular perturbations of the 1-D Laplacian, defined through non self-adjoint 
boundary conditions in the extrema of a bounded interval. For 9q £ C and a, b € R, with a < b 
and b — a = L, the Hamiltonian Hg is defined by 



D(H^ o0 ) = {ueH 2 (U\{a,b}): 



e-^u{b+) = u(b-); e-i 9 °u'{b+) = «'(&-) 



e 



l u(a-) = u{a+)- e-i e °u'(a-) = u'(a + ) J (2.1) 



where u(x + ) and u(x ) respectively denote the right and the left limits of u in x, while the 
notation u' is used for the first derivative. When \6q\ is small, the analysis of H$ follows by a 
direct comparison with Hq (coinciding with the usual Laplacian: — /i 2 Ar). To fix this point, let 
us introduce the intertwining operator Wg defined through the integral kernel 

r +oc l dk 

We (x,y) = J ^(k,x)e-^y (2.2) 
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with ip- (fc, x) denoting the generalized eigenfunctions associated with our model. These are de- 
scribed by the plane wave solutions to the equation 

«,o-fc 2 )V>-(M=0. (2.3) 

For k > 0, one has 

{e*** + R(k)e-^ x , x<a 

A(k)e^ x + B(k)e-^ x , x e (a,b) (2.4) 

T(k)e^ x , x>b. 

Since -ip-(k,x) fulfills the boundary conditions in (|2.ip . the explicit expression of the coefficients 
in (j2~4l are 

with: c+(0 o ) = e^" + e^ , c_(6> ) = e^ 1 - e5 e ° and 



(2.5) 



d(*o,*)=dct C+ f;Kl "J^t-S = C+ (^ 2 e -^- C -(^o)V^. (2.6) 
yc_(6'o)e 1 '- c+(6» )e 4 h D y 

For k < 0, an analogous computation gives 

iT(/c)e l 7r a: , x<a 
A(k)e^ x + B{k)e^ x , x G (a, 6) (2.7) 

with 

f i(fc)=A(-fc), B{k)^e^ a e 2i ^ L B{-k), 

{ (2-8) 
1 f(k) = T(-k), R{k) = e 4l Ti a e 2l Ti L R{~k). 

In what follows we adopt the simplified notation 

4rj u,n_M(*0> k>o ( B(k), k>o 

A{K do) - \ A(k), k < , B{K ° o) ~ \ B(k) k<0, (2 - 9) 

T(k,e )-^ f(kl fc<Qj fl(Mo)-|£ (A)j fc<0 (2.10) 

Lemma 2.1. XTie operator Wg defined by (|2.2I) verifies the expansion 

We -Id = O(\9 Q \) (2.11) 

m operator norm. 

Proof: According to (|2.4[) . (I2.7[) and to the definition (I2.2[) . one can express the integral kernel 
Wg (x, y) as follows 

(A(fc,0 o )-l)e^^— +l (Ol6) (a0 y S(fc,0 o )e- < * (x+3/) ^ 

(T(k,e ) - i)^-v)— + i (bt+oo) ( x ) / (T(k,e ) - i) e **c»-»)_ 



and 



+ l(6,+oo)(x) /" i?(fc,0o)e- 4 * ( ^ ) ^ + l(_ oo , a) (x) / +0 ° R(k,9 )e-^ x+ y^ 
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The previous expression is rewritten in terms of operators: 



W 6o -Id= l^b)^ 1 (A(k, ) + l^PF' 1 B{k, 9 Q )T 

+ l(-oo,a)- F_1 l(-oo,o) (T(k, ) - 1) T + l(b,+oo) •^ r_1 l(o,+oo) (T(k, 9 a) - 1) T 

+ l(6,+oc)-P-7 r_1 l(-oo,0)^(fc, 9q)F + l(-oo,a)P^ ^-(0,+oo)R(k, 0o)F ■ 



dx 



and P denotes 



where T denotes the Fourier transform normalized as Tu(]z) = J R u(x)e ^h) 1 / 2 > 
the parity operator: Pu{x) = u{—x). 

Since the operators P, J 7 , J 7-1 and multiplication by the characteristic function of a set are 
uniformly bounded with respect to 9q, we get: 

\\We -Id\\ <C(\\A(k,0 Q ) - l\\ L oo m + ||B(Mo)IU~(R) + \\T(k,0 Q ) - 1|U»(r) + \\R(k,0 o )\\ L ™ m ) . 

From (|2.8|) . it is enough to estimate for k > the terms at the r.h.s of the inequality above to get 
the L°°(M) bounds. Moreover, from the definition of the coefficients c_(#o), c+{0q) an d d(0o,k), 
we have: 

c_(0 o ) = 0(|0o|), c+(0 o ) = 2 + 0(|0o|), d(9 ,k)=4e- i ¥ +O(\9 \), 

where the upper bound of 0(|#o|) holds with a universal constant. The previous equation gives 
\d(0o,k)\ > 1 when \0q\ is small enough, and using (|2.5p . this implies: 



\A(k,0 )-l\ 



2c+{0 o )e-'— -d(6 ,k) 



\O(\0o\)\ 



\T(k,0 Q yi\ 



\d(0 o ,k)\ 
(c + (9 Q ) 2 -c^(0 o ) 2 )-d(9o,k) 



< C\0 O \ 



\d(0 Ol k)\ 
|2M0 o ) 2 sin(f )| _ \O(\9 \ 2 )\ 



\d(9 ,k)\ \d(9o,k)\ 
\2c.(9 )\ \O(\0 O \)\ 



\d(0 o ,k)\ 



< C\0 Q \ 



\B(k,6 )\ 



<C\6 \, 



\R(k,6 )\ 



\d(9 ,k)\ \d(0 o ,k)\ 
|2 C+ (0 o )c_(0o)siii(f )| _ \(2 + O(\0 Q \))O(\0 o \ 



\d(0o,k)\ 



\d(0 o ,k)\ 



< C\0 Q \ 



□ 

According to the result of Lemma l2.1l for \0q\ small enough, W§ is invertible; in particular one 

has 

We =l + O(\9 \), Wg- 1 = l + O(\0 o \). (2.12) 
Then, it follows from the definition (|2.2[) that Hg are Hq conjugated operators with 

<,o = We Hl Q W e - 1 . (2.13) 

This relation allows to discuss the spectral and the dynamical properties related to Hg for small 
values of the parameter 9q. 

Proposition 2.2. There exists c > such that: for any 0q with \0q\ < c, the following property 
holds: 

1) The operator Hg o has only essential spectrum defined by a ess (Hg g ) = R + . 

2) The semigroup associated with Hg is uniformly bounded in time and the expansion 



e -^K , _ e -itH; 



<i-»+O(\0o\) 



(2.14) 



holds uniformly intE 



Proof: 1) According to (|2.13[) . one has 

«,o -»)~ 1 = We B (H* - zy 1 Wg\ 

with We , Wg' 1 bounded in L 2 (R). This implies: a(H% a fi ) = a ess (H^ ) = K+. 
2) Since Hg o are Hq conjugated by Wg , we have 

e~ ttH %-o Wga = W eo e- itH w. 

For 1 6*o | small, one can use the expansions (|2 . 1 2[) to write 

e -«H? .o = We e- u <°Wi o l = (1 + O(\0 O \)) e" lt ^ (1 + O(\9 \)) . 

Recalling that Hq is self-adjoint (it coincides with the usual Laplacian), and (e~ ltH °<° J is a 
unitary group, one has 



e -HH ao ,c = g-«H£ + O(\6 \). 



□ 



Remark 2.3. It is worthwhile to notice that Hg is not self-adjoint (excepting for 9o — 0) neither 
accretive, since 



Re 



= h 2 Im[[u(a-)u'(a-) - u(b+)u'(b+)] (l - e"^)} , u e D«, ) 



ftas not a /ixerf sign. Thus, it is not possible to use standard arguments to state that e ltH <>o-° is a 



contraction. Nevertheless, for small values of the parameter 8q, the result of Provosition Wl^ allows 
to control the operator norm of e "a- uniformly in time, and states that the time evolution 
generated by Hg g is close to the one associated with the usual Laplacian Hq Q . 



Spectral properties of Hamiltonians obtained as singular perturbations of Hg can be discussed 
using standard results in spectral analysis adapted to this non self-adjoint case. 

Lemma 2.4. Let V = Vi + V2, with V± G L°°((a,b)) and V2 a bounded measure supported in 
U CC (a, 6). Then: a ess {K,fl + V ) = °^( H Lo)- 



Proof: The proof follows from the first point of Corollary 13.41 below in the case 9 = 0. □ 

2.2 Numerical computation of the time propagators for small \9 \ 

This part is devoted to the numerical comparison of the propagators e~ ltH<> °- v and e~ ltH °< v where 
V is a locally supported perturbation of Hg with 

Hl, v =K,o + V, and V € L°°((a, 6)) . 

Using discrete time dependent transparent boundary conditions for the Schrodinger equation, it 
is possible to compute the propagator e~ °> v with a Crank-Nicolson scheme, see 



To compute the propagator e 9 o> v , the key point is to integrate the boundary conditions in 
(|2.1[) in the resolution in a way which preserves the stability. This is performed by integrating the 
boundary conditions in the finite difference discretization of the Laplacian at the points a and b. 
To simplify the presentation, we suppose temporarily that the interface conditions occur only at 
0. So we want to write a modified discretization of the operator with the condition 

f e-*«(0+)=it(Q-) f21g s 
I e-l e °w'(0+) =u'(Q-). V ' ; 
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For a given mesh size Ax, we introduce the discretization of R: Xj — j Ax with j £ Z. For j =^ 0, 
the number Uj will denote the approximation of u(xj), while uq will denote the approximation of 
it(0 _ ) and Uq will denote the approximation of u(0 + ). If the function u is regular on W, we can 
use the usual finite difference approximation 

*\ = Uj -- 1 -2 U , +Uj - +1 
ax- 2 / . Ax^ 

for j ^ {—1,0,1}. Due to the regularity constraint, this approximation is written correctly for 
j = — 1, and respectively for j = 1, only when considering the continuous extension of the function 
from the left, and respectively from the right, which leads to 

d 2 \ u_2 — 2w_i+«o /d 2 \ Uq" — 2tii + "2 

dx^"J_ 1 = Ax 1 ' Vdx^V^ Ax 1 ' 

With the hrst relation in (|2.15[) . the approximation at j = 1 is 

/a \ e 2 m — 2u\ + U2 



- 



\dx 2 J -, Ax 2 



(2.17) 



At j = 0, due to the possible discontinuity of a function u verifying (I2.15[) , we define u~ on R as a 
regular continuation of u^-o^q). More precisely u~ 6 C 2 (R) is such that u~ — u on (— oo,0) and 
we get the following approximation at j = 

—u] = — - 2m ° + (2 18) 

dx 2 / q Ax 2 

This method corresponds to the introduction of a fictive point which allows to write the finite 
difference approximation for 4-^ and to calculate u'(0~) and u'(0 + ) in (|2.I5[) by using the same 
points of the space grid 

I e 2 u ^ = u 

| e~~i e °{ui - u%) =Ui -u . 

The resolution of the system above gives: %" = (1 — e~ e °)uo + e~i B °u\ and (|2.I8p becomes 

U^J = Ax^ • (2 ' 19) 

Therefore, from the boundary conditions in (|2.1j) . the scheme to compute the propagator e~ ltH "o- v 
is obtained by using the modified Laplacian corresponding to the application of (|2.I7|) and (|2.19l) 
at x = a and x = b. When 9o is small, the equations (|2.I7p and (I2.19[) approximate well the 

usual finite difference equation (|2.I6[) . therefore the solution e~ ltHe a> v will be close to the solution 
e ~ ltH o,v ( as expected. 

After the change of variable x' = — 1, where i = the problem for e~ ltHe °' v is the following 

id t u(t,x)= -^d 2 x + V(x) u(t,x), t > 0,x GR\{-1,1} 

e^u(i,-I+) =«(i,-l-), ei e °<9 ;E tt(t,-l + ) =d x u(t,-l-), t>0 ( 2 .20) 

e-T U (t,l+) = u(M~), e-l 9o 9 x u(i,l+) = a x it(t,l-), i>0 

w(0, x) = ui(x), xgR, 

where V^(x) = V((x + l)£ + a) and G C°°(R) is the initial data. The resolution will be performed 
on the bounded interval [—5,5] using homogeneous transparent boundary conditions valid when 
suppu/ CC (—5, 5). 

Set Ax = 4 and consider the uniform grid points Xj = jAx for j G {— 5 J, . . . , 5 J}. Then using 



(|2.17[) and (|2.19[) . the Crank- Nicolson scheme for the modified Hamiltonian is obtained from the 
Crank-Nicolson scheme in 6 25 by replacing the usual discrete Laplacian by the modified discrete 
Laplacian defined below 



I , uj i-(l + e^)«j' + eTiS j+ i, if j = ±J 

Ae ° Mj = 1 e± ^ u j-i - 2u j if j = ±J- 

— 2uj + Uj + i, else . 



The discrete transparent boundary conditions at x — — 5 and x — 5 are those used for e itH o,v i n 

mm- 

For a given time step At and for 9q — * I m $o , we present some comparison of the numerical 
solution u g l o to the system (|2.20p . given at time t n = nAt by the scheme described above, with 
the numerical solution u n to the reference problem, computed by taking 0$ = 0. In particular, the 
numerical parameters are the following: I = 1, h — 0.03, J = 30, At = 0.8, and the comparison is 
realized with the initial condition equal to the wave packet 

t \ ( ( x ~ x o) 2 / X 
ui[x) = exp I — h ik(x - x ) 

where a — 0.2, k — and the center xq will be specified in each simulation. 

Three simulations were performed corresponding to different values of the potential V and of 
the center xq. The first test was realized with V — and xq = —3. Although the comparison 
presented here can be extented to more general potentials, the two other tests were realized in the 
case where V is a non trivial barrier potential 

V = V Q I (a , b) , 

where Vq — 0.8: for this potential one test was realized with an initial condition localized at the 
left of (—1, 1) by taking xq — —3, and the second with an initial condition localized in (— 1, 1) by 
taking xq — 0. The solution Ug g is represented, at different time t, next to the reference solution 
u n in the Figures [TJ [5] and [31 for the fixed small value 6*o = 0.09i. We remark that Ug g has the 
same qualitative behaviour than the reference solution. 

In the case V = 0, the solution Ug o corresponds to an incoming function from the left which 
goes near the domain (—1,1). When time grows, it crosses the interface points and leaves the 
domain. 

In the case of the barrier potential with = —3, the solution is splitted in two parts: a first 
one which passes through the barrier; and a second one which is reflected and goes out of the 
domain. 

In the case of the barrier potential with xq — 0, it appears that the wave packet is splitted in 
two outgoing parts: one which leaves the barrier on the left and the second on the right. The part 
on the right is more important and goes out faster, it is due to the sign of the wave vector k. 

In the three tests described above, although some oscillations occur when crossing the inter- 
faces x = — 1 and i = l, the quantitative comparison gives also good results. In particular, we 
represented in Figure [4] the variation with respect to Im #o of the maximum in time of the L 2 
relative difference: 



\\uS - u II 

D 6a = max 100^ — — - , (2.21 

l<n<N ||u/ M 



where N = 400 is the number of time iterations. It shows that, for every case, the difference tends 
to when Im 8o tends to 0. Moreover, the graphic of Dg is a line which validates the result 
(|2.14[) . We note also that the difference in the case of a barrier potential is smaller then in the 
case V = 0. This may be due to the fact that the error coming from the interface conditions 
is compensated by the exponential decay imposed by the barrier. Then, in the case of a barrier 
potential, we note that the difference is more important when the initial condition is supported in 
(—1,1). It can be explained by the fact that the solution crosses the two interfaces, at x = — 1 and 
x = 1, whereas, when the solution comes from the left, only the interaction with the first interface 
x = —1 is relevant, which is also a consequence of the exponential decay in the barrier. 
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ure 1: Case V = 0: modulus of the functions Ug and u n at different time t with xq = —3 and 
= 0.09i 




ure 2: Case of a barrier potential: modulus of the functions Ug and u n at different time t with 
= -3 and 6 = 0.09i 
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Figure 3: Case of a barrier potential: modulus of the functions Ug and u n at different time t with 
x = and 9 = 0.09i 



Difference with the reference solution 



1.8 




lm(9.) 



Figure 4: Variation with respect to Im 8q of the difference Dg , defined in (I2.21j) . in the case V = 
with xq = —3, and in the case of a barrier potential with xq = — 3 and xq = 
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3 Exterior complex scaling and local perturbations 



Spectral deformations of Schrodinger operators arising from complex dilations form a standard 
tool to study resonances. This technique - originally developed by J.M. Combes and coauthors 
in [3] [TO] for the homogeneous scaling in L 2 (R"): Uei/)(x) — e~2~il>(e x) - allows to relate the 
resonances of the Hamiltonian H = — A + V with the spectral points of a non self-adjoint operator 
H{9) = UeHUg 1 with 9 G C. If the potential V is dilation analytic in the strip {9 G C | |Im 0\ < a }, 
the poles of the meromorphic continuation of the resolvent (H — z) in the second Riemann sheet 
are identified with the eigenvalues of H (9) placed in the cone spanned by the positive real axis 
and the rotated half axis e~ 2llm9 M. + . We refer the reader to [23] for a summary and we recall 
that many variations on this approach have been developed since, see [31] [30] [40] and [29] for a 
short comparison of these methods. In particular for potentials which can be complex deformed 
only outside a compact region, the exterior complex scaling technique appeared first in [S3] in the 
singular version that we reconsider here. Meanwhile regular versions have been used in |31j and 
extended with the so called "black box" formalism in [56j . 

In this section, we consider a particular class of exterior scaling maps, Ug, acting outside a 
compact set in ID and introducing sharp singularities in the domain of the corresponding deformed 
Hamiltonians. Let us introduce the one-parameter family of exterior dilations 

e- e (x-b)+b, x>b 

x, x G (a, b) (3-1) 

e~ e (x — a) + a, x < a . 

For real values of the parameter 9, the related unitary transformation in L 2 (R) is 

e?tp(e e (x - b) + b), x>b 
U e ip(x) = { tp{x), x G (a, b) (3.2) 

e^ip(e e (x — a) + a), x < a . 

Local perturbations of Hg (0) are defined by 

<,v(0)=<,o(0) + V, (3.3) 
with supp V C [a, b]. In what follows we will assume 

V = Vi+V 2 , Vi e L°°((a,b)), V 2 eM b {U) with U CC (a, b) . (3.4) 
Under these assumptions, 

D(KA°)) = i u e # 2 ( R \ {«. b > U }) I i 1 - x)« G D(f$, )0 (0)), ~h 2 u" + V 2 u e L 2 (U) } , 
where x G Cq°((o, b)) and x( x ) — 1 f° r x G U. 

In particular, since V2 is a bounded measure, the domain D(H% o v (0)) is contained in ff 1 (R\ {a, b}). 
The conjugated operator 

H^ v {9) = UoHkyWUi 1 (3-5) 

is defined on D{H^ o V {e)) = |w G L 2 (M) C/^ 1 ?/ G D(^ o V (0)) j. The constraint l^u G D{H^ q V {Qi)) 
compels the boundary conditions 

e-^ 9o+ ^u(b+) = «(&-); e-i( e ° +e )u'(&+) = «'(&-) 

(3.6) 

e -i(e +e) w ( a ^) = e -i(9o+fl) u /( a -) = u '( a +) ; 

to hold for any u G D(Hg y(6*))- Thus one has 

£, (^ ,vW) = { w eiJ 2 (R\{a,fe,[/})niJ 1 (R\{a,6})| . —h 2 u" + V2?i G L 2 (U) } . (3.7) 
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The action of H^ y {9) is 

He o y(0)u = [-h 2 V (x)d 2 x + V] u, n(x) = e -Mi»\(.,»)(x) 



(3.8) 



It is worthwhile to notice that this definition can be extended to complex values of 6. For 9 € 
the Hamiltonian Hg o v {9) identifies with a restriction of the operator Q(9) 



D{Q(9)) = {u e H 2 (R\ {a, b, U}) n H 1 ^ {a, b}) \ -h 2 u" + V 2 u E L 2 (U) } , 
Q{9)u = [-h 2 r](x)d^ +V]u. 



(3.9) 



For particular choices of 9q and 9, the quantum evolution generated by the deformed model H@ v (9) 
is described by contraction maps. To fix this point, let us consider the terms 

Re (u,iHg o v (9)u) ; for u 6 D(Hj} v (9)), an explicit calculation gives 



Re(u,iH& OiV (0)u) 



LH 



{-ih 2 (u(a-)u'(a-) - u{b+)u , {b+)) 
'sm(2Im0) / \u'\ 2 dx. 

jR\<a,b) 



e- 29 - e 



h 2 e~ 2Rc 



i(e+eo) e -f(e+e )^| 
(3.10) 



For 9 = 9q = it, with t £ (0, ^) , the boundary terms disappear, and the r.h.s. of (|3.10p is positive 



Rs(u,tfl£. iV (tt» = /i 2 sin(2r) / |u'| 

'" '" ' JR\(a,b) 



dx > 0. 



(3.11) 



Lemma 3.1. For t G (0, the operator iH^ TV {iT) is the generator of a contraction semigroup. 

Proof: As a consequence of (|3.11[) . iH^ T v (ir) is accretive. Moreover, the propriety a ess [H^ T y («r)) = 
e~ 2lT R + in Corollary 13.41 below, implies i\ e p{H^ T v (ir)) for some Ao > and (iH^ TV (iT) + Ao) 
is surjective. Then, a standard characterization of semigroup generators ([51]. Theorem X.48) leads 
to the result. □ 



3.1 Krein formula and analyticity of the resolvent 

In order to get an expression of the adjoint operator of Hg g v (9), we introduce the following operator 
with two-parameters boundary conditions 



D(Qe u e 2 (d)) = {u€D(Q(6)) 



e -§(0i+%(b+) = e -^+^u'{b+) = u'(b-) 

e -i(»i+fl) u ( -) = u ( +); e-4( fla+M V(o-) = u'(a+) 



Sfli.82 (*)« = [-h 2 ri{x)d 2 x + V] u, r](x) = e~ 2e ^ 



(3.12) 
(3.13) 



where Q{9) is defined in (I3.9[) . Indeed, by direct computation 

i) Qs 1 ,e 2 (9) identifies with the original model Hg o v {9) for the choice of parameters: 9\ = 9o and 
9 2 = 39o 

ii) the adjoint operator (Q,8 1 ,e 2 (@)) is given by 

{Qe^eMY = Q-e 2 ,-eM ( 3 - 15 ) 

Like Hg o v (9), the Hamiltonian Qe 1 ,e 2 (9) is a restriction of the operator Q{9). In this context, we 
fix a boundary value triple {r*? =12 ,C 4 } with P? : D(Q{9)) ->■ C 4 



rfy = h 2 



I -e~IV'(5+) \ 
-Mb-) 

V e-«V(o") 7 



/ e"iV(& + ) \ 

V/(a+) 
V e~iv(a~) / 



(3.16) 
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and (r?,lt) : D(Q(9)*) -> C 4 x C 4 surjective. For all € D(Q(6)), these maps satisfy the 
relation 

(V, Q(0W) L2(m - (Q(6)i>, cp) L2{m = (ifv, r^) c4 - (rf^, rV) c4 (3.17) 

(for the definition of boundary triples and the construction of point interaction potentials in the 
self adjoint case see [27] and [3]). Let A, B e C 4 ' 4 be defined as 



a(0i,0 2 ) 



am, 



2) 



-o(- 




B = 

1 

1 



the boundary conditions in (|3.12p arc equivalent to 

Ar^V = BY 6 2 ij) . 

Let H^ D v {9) denote the restriction of Q(0) corresponding to the boundary conditions: Tii[> 
0; this operator is explicitly given by 



(3.18) 
(3.19) 

(3.20) 



2 \D 



—h A 



(a,b) 



V 



-h 2 e- 2e A N 



(b,+oo)- 



Its spectrum is characterized as follows 

° {H%d, V {Q)) = e- 2e M+ U a (- 



^ 2 Af Q , b) + V 



(3.21) 



(3.22) 



It is possible to write (Qe 1 ,e 2 (0) — z) 1 as the sum of (H^ D v (9) — z) plus finite rank terms. 
Such a representation will be further used to develop the spectral analysis of Hg v (9) where 
our Krein-like formula will allow explicit resolvent estimates near the resonances. The space 
Nz,0 = Ker(Q(0) — z) is generated by the linear closure of the system {ui iZ \- x where u^ z are the 
independent solutions to (Q(0) — z)u = 0. The exterior solutions to this problem, itj jZ , i = 1,4, 
are explicitly given by 



u\ >x {x) 



-(x-b) 



Ui, z {x) = l(_ 00:a )e 



- (x—a) 



(3.23) 



where the square root branch cut is fixed with Im^A > 0. This assumption implies Im V ze 20 > 
for all z e C\e- 2tIm9 R + . The interior solutions, Ui lZ , i = 2,3, can be defined through the following 
boundary value problems 

f [-h 2 d 2 x +V-z) u 2<z = 0, 

1 u 2 , z (a) = 0, u 2 , z (b) = 1 



[-h 2 dl +V-z] u 3 ,2 = 0, in (a, b) 
u 3 , z (a) = 1, u 3 . z (b)=0, 



(3.24) 

with z G C\ap(H^ D v (6)). Owing to the property of the interior Dirichlet realization in (a, b), the 
are unique and locally H 2 near the boundary. We consider the maps: 7(-, z, 6) = 



solutions 



-1 



(a,b) 



(■^i|jV e ) ' W ^ n ^iIaT denoting the restriction of T® onto Af z ,$, and q(z,9,V) — T 2 > 'y(-, z,9). 
These form holomorphic families of linear operators for z in a cut plane C\e- 2lIm9 R+. Their matrix 
form w.r.t. the standard basis {&,■} . _ 1 of C 4 and the system {ui^ z } i=1 is: -fij(z,9) = Ci(z,9)5ij 
with 



ci{z,9) = 



1 



1 



hV 'ze 26 



> c 2 



and 



/ ihe 



q{z,6,V) = - f? 



-ft?' C3 = V' C4(z ' 



-u' 2 ,M u' 3 Jb) 
-u' 2 Ja) u' 3 (a) 



le 2 



hV ze 29 



ihe 



(3.25) 



(3.26) 
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Lemma 3.2. Let p £ L 2 (R) and j — 1, ...,4; the relation 



A{H h NDy (e)-z) \ =( 1 (e J ,z,0),p) 1 



holds with 9 £ C and z £ p {H^ D v (0)) ■ 

Proof: Let: / = (H'^ DV (9) - zy 1 ip. This function is in D{H% D V (6)) so that: 
Q(9)f = H% DV {6)f and T e J = 0. The l.h.s of (13~2T1) writes as 



(3.27) 



T 9 2 {H h NDy {9)~z) \}=(e j ,Tlf). 



Since e ? - = Ff j{§.j, z, 9), we have 

= (t(%,*,0), QW/) ia(R) - (Q(9)'y^z,S), /) i2(R) 

By definition, 7(ty-,5, 0) 6 A4 and the r.h.s. writes as 

(7(^,2,0), Q(8)f) L z W ~ (Q(0h(£j,z>0), f) LHm = He^O), (H h NDy {9) - z) /) 

= (7(^,2,0), ^) i2(R) . 



L 2 ( 



□ 



Proposition 3.3. The resolvent (Qe 1} e 2 (^) ~ z ) 1 allows the representation 

{QeuoM-zT 1 = {H h ND>v {0) - z y l 

4 



(3.28) 



and one raas: a- e5S (fl<?i,0 2 (0)) = cr ess {H% Dy {6)) = e W R + . 

Proof: Let us consider the r.h.s. of this formula: the operator (H^ D v (9) — z) 1 is well defined for 
z e C\a(H^ Dy (9)). The vectors 7(e i; z,6»), i = 1,...,4, are given by (j3~2lj) . (|3T23]l . (j3~23|) . while the 
boundary values of 2 - appearing in the definition of the matrix p. 261) - are well defined whenever 
z e C\o- p (H$f Dy (9)). Therefore, the r.h.s. of (1335)) makes sense for z £ C\ {o-{H% Dy {9)) U %} 
where 7o is the (at most) discrete set, described by the transcendental equation 



det (Bq(z, 9, V) - A(0i, 2 )) = 0. 



(3.29) 



It is worthwhile to notice that C\ {a(H^ D v {9)) U 7o} is not empty. Let us introduce the map 
R z (ip) defined for <p £ L 2 (R) by 

Rz((fi) =<j>-ip, 
<P={H h NDy {9)- Z y 1 p, 

4 

yj= J2 [(B^,^V)-A)- 1 s]^(7(e J ,z,^),^) i2(R) 7(e l ^^), 

i,J = l 

with q(z, 9, V) given in (I3.26P and z £ p(H'^ D y (0))\7o- In what follows we show that: 

Rxifp) = (Qe lt eM - z)~V Since H^ D V {9) C Q(0) and 7(^,2,0) e Ker(Q(0) - z), one has: 



(4d,v(8) - z ) V 7fe, 2, 9) £ D(Q(9)) . 
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This implies R z (f) G D(Q(9)). To simplify the presentation, we will temporarily use the notation 
q = q(z, 8, V). Being <fi £ D(H^ D v (8)), we have: T\<j) = and the following relation holds 

(Bq - A) r? {4> - ,/,) = - (BT e 2l ( :i z, 8) - A) = (-BT* + AT?) i>, (3.30) 
where tp € A4 e and 

7(-,M)r?k.,. = i 

have been used. At the same time, the ro-th component of the vector at the l.h.s. can be expressed 
as 

[(Bq - A) r? (<f> - ^)] „ = [- (Bq A) rfy] n 



J2 [(Bq-K)T e ll (e l ,z,e)} n \(Bq-A) 1 b] (^,2,0),?) 



L 2 (R) 



Recalling that Tfj^, z, 8) = e», we get 
[(Bg-A)r? (0-^)]„ = 

4 

- £ (Bq~A) rn \(Bq — A) -1 sl ( 7 ( % ,z,0),^) 



J = l 



Taking into account the result of the Lemma 13.21 this relation writes as 

(Bq — A) r? (cj>-T])) = -BT e 2 <j>. 



(3.31) 



Combining (j3~30l) and (j3~3Tj) . one has: -ATfip = BY e 2 (<f> - tp), and, adding the null term ATfcp at 
the l.h.s., 

ATlRzitp) =BT e 2 R z (<p), 

which, according to (|3.20[) . is the boundary condition characterizing Qg lt g 2 (8) as a restriction of 
Q(8). Then we have: R z (<p) & D (Qg u g 2 (8)). Furthermore, 



(Qe u e 2 (d) - z) R z (p) = (Q(0) - z) R Z ( V ) = ip- (Q(8) -z)rP = tp, 



(3.32) 



where (Q(8) — z)^y(e i ,z, 8) = has been used. This leads to the surjectivity of the operator 
(Qe 1 ,0 2 {Q) ~ z ) f° r an Y {<t(H'^ d v (8)) U 7o}- The injectivity is obtained using the adjoint 

of (Qg lt g 2 (9) — z )- Indeed, the equality (|3. 151) implies 

which, from the result above, appears to be surjective for all z such that z G C\ ^(H'^ D v (8)) U 7o|, 

where 7o is the discrete set of the solutions to (|3.29p when replacing: 9 = 9, 9\ = —8 2 and 8 2 = —9\. 
As a consequence of (|3.22l) . we have 



* e o- (i?& AV (0)) ^ a (H h NDy (8)) . 

It follows that for any z e C\ {cr^^.v^)) UT}, where T = 7B U | 
(Qei.e 2 (^) — z ) i s surjective and 



z s.t. z € 



To}, 



the operator 



Kei(Q ei ,eM -») = [Ran (Q ei ,e a (0) - z)*) X = {0}. 
Wet get that Vz e C\ {o-(iJ^ DV (l?)) U T}, (Qe 1 ,e 2 (°) ~ z ) is "avertible and (Q ei ,e 2 (0) - z) _1 = 



Moreover, for such a complex z, the difference R z — (H^ D v (8) — z) is compact. Then, we 
conclude that 
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(for this point, we refer to [52], Sec. XIII. 4, Lemma 3 and the strong spectral mapping theorem), 
and the equality (I3.28P holds as an identity of meromorphic functions on C \ e _2e R + . □ 
As a direct consequence of the previous proposition and the identification (I3.14[) . the represen- 



tation of the resolvent ( Hg a v {8) — zj is obtained by replacing the matrix A in (|3.28[) by the 
matrix A — A(8q, 38q) where the matrix A(9i, 82) is given in p,18[) . It follows that for the matrices 



■e ? 



A = 



h 2 



-e 2 



e 2 



V 



e 2 



V 



B = 



1 

J 



(3.33) 



the result below holds. 



Corollary 3.4. Let V and A, B be defined as in j3.4\ ) and p. 33 

allows the representation 

{Hty^-z)- 1 = (H* Dy {6)-zy 1 

4 

T \{Bq{z,6,V)-A)- l B\ ( 7 ( %I z, 0), •) 

L J ij 

i 

and one has: a ess 



The resolvent I Hg v (8) — z 



L' 2 ( 



(3.34) 



<(Hly{e))= rr. 



-20T1 



The analyticity of the resolvent (^Hg g v (8) — z^j w.r.t. 8 is an important point in the theory 

of resonances. The former is obtained in the next proposition as a consequence of the formula 
(I3.34p . the latter is developed in the next section. 



Proposition 3.5. Consider Oq G C, V fulfilling the conditions \3.4\j , and let S a = {9 G C | |Im 8\ < a } 

for positive a < ?• Then, there exists an open subset O C {z G C |Re z < 0} such that V# G S a , 

O C p(Hl v (6)). Moreover, Vz G O, the map: 8 » (h'^8) - z)" is a bounded operator- 
valued analytic map on the strip S a . 

Before starting the proof, we note that this result implies that the 8 dependent family of 
operators Hg g v (6) is analytic in the sense of Kato in the strip S a (definition in (52)). 
Proof: The equality (|3.34[) holding as an identity of meromorphic functions on C \ e 



- 2e K+, the 

poles on the l.h.s. identifies with those on the r.h.s. . For 8 G S a , the characterization (|3.22l) 

zV> 



implies that the map z h-> (H^ D v (8) — z) "is analytic when Re z < 0. It results 
°~ ( H e ,v( 6 )) n {z e C |Rc z < 0} = % n {z g C |Re z < 0} , 

where 

T = {zeC |det (Bq(z, 8,V) - A) = 0} . 

\fze 6 , and therefore q(z,6,V) = q(z, 0, V), 



,28 



Noting that for Re z < and 8 G S a we have vze' 
we get: V6> G S a 

a(H e l oy (8)) n {z e C I Re z < 0} = {z G C | Re z < and det (Bq(z, 0, V) - A) = 0} , 

which is a discrete set independent of 8. This implies that there exists an open subset O C 
{z G C I Re z < 0} such that V0 G 5 Q , O C p f#£ oV (0)) . 

In what follows, we fix z G O. The equation (|3.34p gives, for any 8 G S a 
-1 ,^ \-i 



«,v(f)-*) 



{ H NDy{8) 



[(Bq(z, 0, V) - A)" 1 fl] . ( 7 (e r z, 0), -) i2(R) 7 (e i; z, ( 



(3.35) 
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where we used q(z, 9, V) = q(z, 0, V), and we want to study the analyticity of the r.h.s. with respect 
to 9. 

Let us start with the operator H^ D v (9): V(9 G S a , it is a closed operator with non empty 
resolvent set. Moreover, D (H% Dy {6)) docs not depend on 9 and Vip £ D {H% Dy {9)) , V/ G L 2 (K) 
the map 

e^{f,H h NDy m) L , m 

is analytic on S a . This means that H^ D v (9) is analytic of type (A) following the definition in 
|52) . In addition, since Re z < when z G O, (j3~22j) implies z G p{H% Dy {6)) for all (9 G S a . 

Then, it results from the analyticity of type (A) propriety that the map 9 M> (H^ D v (9) — z) 1 is 
analytic on S a . 

Concerning the finite rank part in (|3.35p . for any z with Re z < 0, the functions j(ei, z,0), 
i = 2, 3, given by (|3.24p(|3.25p . do no depend on 9, and the functions 

- ii - ii 

%e 2 ■ V ze 2t) (x-i)) %e 2 ■ Vzci 2S (x-a) 

7(ei,z,g) = ; ^_ l (b , +oo) e' " , i{e A ,z,9)= - 7 == l ( _ 0o , a) e 1 



2P 



are such that V/ G L 2 (R), 9 n- (/, 7(e i; z, 9)) L2 ^ is analytic on S a . It follows that \-> 7(6,, z, 0) 
is a L 2 (R)-valued analytic function on S a . This propriety holding for any z with Re z < 0, we have 
also that 9 i-> 7(ej,z, #) is a L 2 (R)-valued analytic function on 5 Q . Therefore, for i, j = 1, ...,4, 
the operator with kernel 7(^,2;, 0) (g> 7(e^,f, 0) is analytic w.r.t. on the strip S a . It allows to 

conclude that 9 t— > ^Hq y{9) — z ) ^ s a bounded operator- valued analytic map on the strip S a . 

□ 

3.2 Resonances 

Next consider: Hg g v = H' e l g + V, with V fulfilling the assumptions (|3.4p . Local perturbations of 
Hg Q q can generate resonance poles for the associated resolvent operator. These can be detected 
through the deformation technique by means of an exterior complex scaling of the type introduced 
in (|3.2p . To fix this point, let us introduce the set of functions 

A = \u u(x) =p{x)e~ fSx2 , (3 > 0} , (3.36) 

where x G M and p is any polynomial. The action of Ug on the elements of A is 

' e e2 p{e e {x-b) + b)e- fi ( e ' >{x ~ b)+b )\ x>b 
U e u(x) = { p(x)e- 0x \ x€(a,b) (3.37) 

e%p(e e {x - o) + a) e -0(«*(»-«)+<») a , x < a . 

If Re(e 2e a; 2 ) > ex 2 for some e > 0, the function Ug u belongs to L 2 (R). In particular, for all positive 
a < j, the map 9 H> Ug u is a L 2 -valued analytic map on the strip S a = {9 G C | \lm0\ < a}. 
According to the presentation of [35] , the quantum resonances of Hq q v are the poles of the mero- 
morphic continuations of the function 

z i-> F(z,0) = (u, (Hj} o V ~ z)~ 1 v) , u,veA, (3.38) 
from {z G C, Imz > 0} to {z G C, Imz < 0}. 

Proposition 3.6. Let 9o G C and u,v G .4. Consider the map: z 1— > i^(z,0), a.e. defined in 
{C I Imz > 0}. 

XTie map z 1— ^ _F(z,0) ftas a meromorphic extension into the sector {argz G (— J, 0)} o/i/ie 
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second Riemann sheet. 

2) The poles of the continuation of F(z,0) into the cone 

K T = {argz e (— 2r, 0)} , with r < —, 
are eigenvalues of the operators H@ o v (9) with r < Im 9 < j. 

Proof: 1) The proof adapts the ideas underlying the complex scaling method (see [23]) to the 
particular class exterior scaling maps Ug. 

Consider the strip S a = {9 6 C | |Im#| < a} for a positive a < j. Then, consider the corre- 
sponding set OcC given in Proposition ^. 51 and fix z G 0. According to Proposition 13.51 and to 
the properties of Ugu, u € A, the function 



F(z,e) = (u s u,(Hl >v (9)-z) 'Ugv) 



is analytic in the variable 9 G S a . When the exterior scaling Ug is an unitary map and one 

has 

F(z,6) = F(z,Q), V6>GS Q nR. 

Since F(z,9) is holomorphic in 9 and constant on the real line, this is a constant function in the 
whole strip S a and 

F(z,6)=F(z,0), V9eS a . 

Now, fix 9 G S a such that Im 9 > 0. It follows from Corollary I3H that a ess (H e h o V (9)) = 
e -2iim e R+ and the map z ^ F(z,0) is meromorphic on C \ e ~ 2ilm e R+ such that 

F(z,9) = F(z,0), VzgC (3.39) 

Since z i-> 0) is meromorphic on C\K+, the equality p.39p holds as an identity of meromorphic 
functions Vz with Im z > 0. We conclude that F(z,9) defines a meromorphic extension of F(z,0) 
from the set {z G C |Im z > 0} to the sector {argz G (— 21m 9, 0)}. 

2) Consider < r < f. From the previous point, when 6* varies in r < Im6> < j, the 
maps z i — ^ -F(z, 0) coincide in i^T- with meromorphic extensions of -F(z, 0). Therefore, the poles of 
z i— > F{z, 9) in K T do not depend on 9 and correspond to the poles of the meromorphic extension of 
F(z, 0) in K T . The vectors Ugu, u G A, being dense in L 2 (M), the poles of z n- F(z, 9) corresponds 
to eigenvalues of Hg v {9). □ 

3.3 A time dependent model 

Consider the non-autonomous model Hg o yn\(9o), where V(t) is a family of self-adjoint potentials 
composed by V = Vi + V2 

n 

Vi{t) GC (R + ;i oo ((a,6))), V 3 (t) = - c,) , (3.40) 

i=i 

with {cj} C (a,b), ctj(t) G C 1 (M+;]R). According to the specific feature of point perturbations, 
the domain's definition, given in (13.7[) . can be rephrased as 



D(H^ o V{t) (9 Q )) = {ue H 2 (R\ {a, b, Cl , c„}) n H\S\ {a, b}) \ 

h 2 [u'(cj)-u'(cj)] = ajfyuicj) and (JSU) holds} , (3.41) 



where time dependent boundary conditions appear in the interaction points Cj. 
Most of the techniques employed in the analysis the Cauchy problem 

{idtu = H{t)u 
u t=0 = u 



(3.42) 
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for non-autonomous Hamiltonians, H(t), require, as condition, that the operator's domain is inde- 
pendent of the time (we mainly refer to the Yoshida's and Kato's results ( [36] j37] [S] ) ; an extensive 
presentation of the subject can be found given in In the particular case of Hg Q y( t )(#o); one 

can explicitly construct a family of unitary maps Vt,t such that 

VttoBl >v(f) {9 )V t -£ (3.43) 

has a constant domain. To hx this point, let us introduce a time dependent real vector field g{x, t) 
and assume 

n)g( Cj ,t)=0, j = l, n ,Vt, (3.44) 
hi) supp#(-,t) = U" =1 I C3 ., with: I Cj = (cj - ej,Cj + tj) and Dj I C] = . 

For ej small, g(-,t) has support strictly included in (a, b) and localized around the interaction 
points x = Cj. According to i), g(-,t) satisfy a global Lipschitz condition uniformly in time: then 
the dynamical system 

f Vt = g(vt,t) 

(3.45) 

I Vto = x 

admits a unique global solution continuously depending on time and Cauchy data {to,x}. Using 
the notation: y t = F(t,to,x), one has 

\F(t,t Q ,x 1 )-F(t,to,x 2 )\<e M \ t - to \\x 1 -x 2 \, with: M = sup d y g(y,t). (3.46) 

XEM., t>0 

Consider the variations of F(t,to,x) w.r.t. x: z t — d x F(t, to, x). From (|3.45p . one has 

f i t = dig(F(t,t ,x),t)z t 

< (3.47) 

I z to = 1 , 

with dig denoting the derivative w.r.t. the hrst variable. The solution to this problem is positive 
and explicitly given by 

d x F{t, t , x) = A *9W:to,x),.)d. >Q \/x e M . (3.48) 

According to ((3T48J) , F(t,to,-) is a C°°-diffeomorphism and the map x i— > F(t,to,x) is a time- 
dependent local dilation around the points Cj. In particular, it follows from the assumptions i) and 
ii) that 

F(t,t ,Cj) — Cj, F(t,to,x) — x for all x € R\suppg, (3.49) 

and 

F{t, t , F(t , i, x)) = F(t, i, x) . (3.50) 
The unitary transformation associated with the change of variable x — > F(t,to,x) is 

f (V to , t u)(x) = {d x F{t,to,x)) i u(F(t,to,x)) 

< (3.51) 
Regarded as a function of time, t H> Vt 0i t is a strongly continuous differentiable map and one has 



dtV to ,t = Vto,t 



- (d y g) + gd. y 



2 y-^vii) t y^v 

The form of the conjugated Hamiltonian Vt.t Hg o yu\{6o)Vt ,t follows by direct computation 

Th (a WT — 1,2 /„A \p, i2q , fa .\ „21 



(3.52) 



Vt,t Hl,V(t)(Wo,t = -h 2 v(y) [d y b 2 d y + (d y ab) - a 2 



Vi(F(i , t,y),t)+J2 <*j (t) b( Cj ,t) S(x - Cj ), (3.53) 
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b(y,t) = etio d ^ s < t ^ ds , a(y,t) = hHy,t) [ d 2 1 g(F(s 1 t,y),s)b(y, s )ds . (3.54) 

1 Jta 

After conjugation, the boundary conditions in the operator's domain change as 

h 2 b 2 ( Cj ,t) [u'(cf) - u'(cj)] = b( Cj ,t) a 3 {t) u( Cj ) . (3.55) 

Assume, for all j = l...n 

iv)aj(0)^0; v) aj(0)aj(t) > 0; vi) aj(t) G C 1 (M+;M) . (3.56) 
One can determine (infinitely many) g(-,t) € C° (R+;Co°(R)) such that 

d 1 g(c j ,t) = %@-. (3.57) 

This implies 

QLj (0) 

and (|3.55p can be written in the time-independent form 

h 2 [u'(cj) - u'(cj)] = aj (0) u( Cj ) . (3.59) 
Thus, the Hamiltonians Vt.t Hg g v ^ (#0)^0,* have common domain given by 

Y = {«£ ff 2 (M\{a,6,c i })nF 1 (M\{a,6})| , [335]}. (3.60) 
Consider the time evolution problem for Hg o v ^(9a) 

( id t u = H^ m (9 )u 

{ U ta = U . 



Setting v = Vt t t u, and 

"1 



A(t) = Vt,t o Hl yit) (0 o )V to , t -t 
with D(A(t)) = Y, one has: v to = u t(> , 

( dtV = —iA(t)v 



2 ( d v9) + f}dy 



V to = u 



(3.61) 



(3.62) 



(3.63) 



Proposition 3.7. Let V(t) = Vi(t) + V2(t) be defined with the conditions {3.40fy and 13.56]) . 

Assume: 9 = 9q = ir, with r £ (0, There exists an unique family of operators Ut, s , < s < t, 
with the following properties: 

a) Ut lS is strongly continuous in L 2 (M.) w.r.t. the variables s and t and fulfills the conditions: 
U SjS — Id, Ut tS o U Str — Ut.r for r < s < t and \\Ut, s \\ < 1 for any s and t, s < t. 

b) For u s G D(Hg o y( s )($o)); has: Ut. s u s G D(Hg o v(t)(^°)) f or a ^ ^ — s ^ n particular, for fixed 
t, Ut, s is strongly continuous w.r. t. s in the norm of ff 2 (R\ {a, b, cj}) n iJ 1 (R\ {a, b}). While, for 
fixed s, Ut, s is strongly continuous w.r.t. t in the same norm, except possibly countable values of 
t > s. 

c) For fixed s and u s £ D(Hg Q y( s )(#o))> the derivative ^Ut )S u s exists and is strongly continuous 
in L 2 (R) except, possibly, countable values t > s. With similar exceptions, one has: 4zUt }S u s = 
- iH 6 ,V(t)( e oWt, s U s . 

d) Additionally, ifVx{t) G C 1 (R+, L°°((a, b))) n C° (R+, W^°°(R)) , and aj (t) G C 2 (R+;R), then 
the conclusion of point (c) holds for all t > s without exceptions. 
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Proof: Since the Cauchy problems (|3.6ip and (|3.63[) are related by the time-differentiable map 
Vt,t , it is enough to prove the result in the case of A{t). Let assume the conditions p.44[) and 
(I3.58|) to hold, and start to consider the properties of this operator. 



I) As already noticed (see relation (|3.11|| ). iHg 



flo.V(i) 



is an accretive operator. This property 



extends to Vt r t f^e v(t)(^°)) ^*o,i> which is unitarily equivalent to an accretive operator, and 
to A(t), since, as a straightforward computation shows, the contribution i [| (d y g) + gd y ] is self- 
adjoint. The spectral profile of A(t) essentially follows from the properties of Hg o V (t)(^o)- Indeed, 



we notice that: 



^*>*o-^e ,v(t) 



(0o)V t 



a (Hq v({)(^o)J i since the two operators are unitarily 



equivalent. Moreover, the term i [| (d y g) + gd y ] is relatively compact w.r.t. Vt t t Hg o v ^{do)Vt .t, 
since it has a lower differential order (see definition (I3.62p ). Then, <r ess (A(t)) — a ess (Hg o v ^ (8 j) = 
as it follows from Corollary 13.41 and (A(t) — z) 1 is a meromorphic function of z in 



-2i Im 9q J 
-2ilm 6»n 



C\e- 28lme °M + . This result yields: a ess {iA(t)) = e l{ - f - 21me o)R + aiK } : p(iA(t)) As a 

consequence of the above, one has: i) iA(t) is accretive; 2) — Ao G p(iA(tj) for some A > 0. Then, 
for any fixed t, iA(t) is the generator of a contraction semigroup, e~ lsA ^ , on L 2 (R) (see [ST], Th. 
X.48). In particular, the operator's domain D(iA(t)) — Y, defined in (|3.60l) . is invariant by the 
action of e~ lsA ^ =>■ Y is iA(t) — admissible for any t. Moreover, as iA(t) is the generator of a 
contraction semigroup, from the Hille-Yoshida's theorem it follows that 



I- C p(iA(t)) 
(»A(t) + A) _1 



< i for all A > . 



In particular, for any finite collection of values < t± < ... < i/., one has 



n-=i {iA^ + xy 1 



< 



1 

A* ' 



(3.64) 



(3.65) 



This relation implies that iA(t) is stable (with coefficients M = 1 and (3 = 0, according to the 
definition given in |37]). 

II) Consider the action of A(t) on its domain. Here, the Hilbert space Y is provided with the norm 
H 2 (M\ {a, b, ci, c„}) n iJ x (R\ {a, 6}). For it e F, we use the decomposition: u = ls,\( a .b) u + 
l(a,6) u — u ext + ""i rl . Recalling that the functions g(-,t), b(-,t), a(-,t) are supported inside (a, 6), 
from the definition p.53[) . (|3.62p . one has 



with 



A{t)u in = 



+ 



(A(t) - A(s)) u = (A(t) - A(s)) m 



-h 2 d y b 2 d y + J2 Kcj,t) S(x - 

3 = 1 



-h 2 ((d y ab) - a 2 ) + Vi(F(t 0) t, y), t) - i - (d v g) + gd, 



(3.66) 



(3.67) 



Taking into account the boundary conditions (|3.59[) , the second order term in Q3.67P writes as 

n 

-h 2 (dyb 2 d y ) u in = -h 2 (d y b 2 ) d y u m - h 2 b 2 A^^.yUin -}]aj(0) b 2 (cj,t)u in (cj) S(x ~ cj) , 

3=1 

where, according to ()3.58[) . 



^2 <*j (0) b 2 (cj , t) u in (cj )S (x - Cj) = ^2 otj (t) b{cj , t)u in (c,- ) S(x - Cj) 

3 = 1 3=1 
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From these relations, it follows 

A(t)u in = -h 2 b 2 A (a:b )\ {Cj} u m - h 2 [(dyb 2 ) d y + {d y ab) - a 2 

'1 



+ V 1 (F(t ,t,y),t)u i . 



2 ( d v9) + gdy 



Therefore, we have 

(A(t) - A(s)) u = -h 2 [(b 2 - b 2 s ) A ( a, 6 )\ { c 3 } + (dib 2 - d!b 2 ) d y + (d iat b t - d ias b s ) - (a 2 - <£)] 

1 



+ (Vi(F(t ,t, y),t)- Vi(F(i , s, y),s)) m 



(dig t - dig s ) + (g t - g s )d y 



and 



\\(A(t) - A(s)) u|| i2 ( R ) <C h M(t,s) IMI ff 2 (K \ {a b Cj } )rW i (R \ {Q b}) , 
where Ct is a positive constant depending on h, while 

M(t,s) = \\Vi(t) - Vi(s)|| Loo:((0ib)) + \\a 2 -a 2 s \\ c0{[aM) + \\dia t b t - 5ia s 6 s || ioo((a b)) 



(3.68) 



£ (11^?-^; 



1=0,1 



i 1,211 

»llL«((o,6)) 



lL=o((o,6)) 



From our assumptions, t H> g t ,6 f ,a t are continuous C°°-valued maps, and 1 1-> Vi € C°(R + ,L°°); 
this yields: lim s _>t M(t, s) = 0, and, due to (|3.68p . 



lim||(^(t)-A(«))«|| £afR) =0, 



This gives the continuity of t h-> j4(t) in the £(y, L 2 (R))-operator norm. 

Ill) Consider the map: S'(i) = (zA(t) + Ao) for Ao > 0. According to (|3.64[) . S(t) defines a family of 
isomorphisms of Y to L 2 (R). Moreover, making use of the relations (|3.53[) . (|3 .62(1 and the condition 
13.58[) , the time derivative of S(t) writes as 



-S(t) = -h 2 [d y (2bd t b) d y + {d y , t ab) - 2ad t a] + £ 



" 2a i (t)o/ i (t) 



3=1 



"3(0) 



Cj) 



+ 9tVi(F(t , t, y),t) + (d t F(t , t, y)) d{Vi (F(t , t, y), t) - i 



2 ( 9 v,t9) + (dtg) d y 



where the term dtF(to, t, y) can be written in the form 

d t F(t ,t, y) = -b^iy, t) g{F(t ,t, y),t), 

as it follows by using: F(t, io, F{to, t, y)) = y and the definition of b(y,t). Since the variations of 
the functions a, b,V±, g w.r.t. both the variables t and y arc supported on (a, 6), each contribution 
to iS(t) acts only inside this interval. Therefore, u G D (^jS'(t)) is not expected to fulfill any 
interface condition at x = a, 6, while in the interaction points Cj, one has 



h 2 2b( Cj ,t)d t b(c J: t) [u'( c t) - u'(cj)] 



2a,(*K(t) 



which, according to (|3.58p . reduces to (|3.59p . It follows that: y C D (4rS(t)), and one can 
consider the action of -iS(t) on Y. Proceeding as in point II and using the stronger assumptions: 
Vi(t) G C 1 (R+,L°°((a,&))) n C* (R+,IF 1 ' 00 (IR)), tti (t) G C 2 (R+;R), one shows that f t S(t) is 
strongly continuous from Y to L 2 (M.) . 

Finally, the points I and II resume as follows: i) the Hilbert space Y is A(£)-admissible for all t, 
ii) A(t) define a stable family operators ^-continuous in the C(Y, L 2 (R))-operator norm. Thus, the 
Theorem 5.2 in |37) applies and provides a strongly continuous dynamical system Ut,t i < t-, for 
the Cauchy problem (I3.63[) with ||Wt,t || < 1 and such that: 
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1) For fixed to, the family Lit,t defines an a.e. strongly continuous flow in Y; for fixed t, Ut.t is 

strongly continuous w.r.t. to in Y. 

2) For y EY, the derivative -^Ut,t a y is a.e. strongly continuous in L 2 (R) and one has: 4Mt,t u s = 

—iA(t)Ut,t u s . The conclusions (a), (b) and (c) of the statement follows by using the equiv- 
alence of the problems (13.61)) and (|3.63l) through the maps Vt,t ■ 

Moreover, assuming 

Vi(t) G C* 1 (R+,L°°((a,&))) nC° (R+, W^fR)) , ay(i) G C 2 (R+;R) 
it follows from point III that: S(t) is a family of isomorphisms from Y to L 2 (R), with 

S(t)A(t)S- 1 (t) = A(t) 

and such that S(t) is strongly differentiable. In this case, Yoshida's Theorem applies (see Theorem 
6.1 and Remark 6.2 in [37]) and the conclusion (<£) of the statement follows. □ 

4 Exponential decay estimates 

With this section, with start the analysis of the parameter dependent quantities as h 0. The 
exponential decay is specified with a good control of the prefactors which behave like pj-. These 
estimates are written for potentials with limited regularity assumptions in order to hold for the 
modelling of quantum wells in a semi-classical island with non-linear effect. Some preliminary 
estimates are reviewed in Appendix [X] 

4.1 Exponential decay for the Dirichlet problem 

Consider V G L°°((a, b); R), and W h = + W% an /i-dependent real-valued potential with: 

d [a , b] <V, \\V\\ L °o<-, \\Wt\\ L oo<-, \\w£\\ Mt <-h, (4.1) 
c c c 

where denotes the total variation of the measure /« G A^b((a, b)). The constant c denotes 

a fixed positive value that can be chosen small when it is required by the analysis. We suppose 
that Wi and W 7 ^ 1 are supported in the domain Ut — {x G (a, b); d(x, U) < h} where U is a fixed 
compact subset of (a, b). 

After introducing the differential operators on (a, 6) 

P h :=-h 2 A + V and P' 1 := —h 2 A + V — W h , 

two Dirichlet Hamiltonians are considered 

H h D := -h 2 A + V = P h , with = iJ 2 ((a,6))ni/ 1 ((a,6)), (4.2) 

iJ^, := —h 2 A + V — W h = P h , with D(fl^) = {u G H^((a,b)), P h u e L 2 ((a,b))} .(4.3) 

For a real energy A G R we consider the Agmon degenerate distance associated with V 



d Ag (x,y,V,X)= y/(V(t)-X)+ dt, x<y. 

J X 

And an other tool that will be useful here is the /i-dependent H k norm 
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Proposition 4.1. i) Consider / = /1+/2 with fx £ L 2 ((a,b)) and $2 £ Mb{{a,b)). IfV—Re z> c 
with \\V\\l^ < -, then any solution u £ Hq((q, b)) to (P h — z)u = f satisfies 

h 1 ' 2 sup \e*&u(x)\ + \\he%v!\\ L * + \\e*u\\v> < ^ ( \\A\\l* + ^tIMm,) , (4.5) 

xe[a,b] "V « 2 / 

with (f(x) = dA g (x, K, V, Re z) and K D supp /i U supp fi . 

ii) Consider f £ L 2 ((a,b)). If V — Re z > c with ||V||£°c < then any solution u £ D(Hp) to 
(Hp — z)u = f satisfies 

h 1 ' 2 sup \e^u{x)\ + WheWWv + \\e%u\\ L i < (||u|| £ . + , (4.6) 

x€[a,b] " 

with <p(x) = d(x, K' U U, V, Re z) and K' D supp / . Especially when z ^ a{H^j), we have: 

h 1 ' 2 sup \e^u{x)\ + \\he*u'\\ L 2 + \\eJu\\ L , < %^ ( ) + l) , 
xe[a,b] h \d{z,a{H^)) J 

and, when z = E h is an eigenvalue of H^, the related normalized eigenvector satisfies 

1/9 | / xi 11,^/11 ii V it Ca b c 

h < sup |e <> u(x)| + \\he h u \\ L i + \\e h u\\ L 2 < — ^- , 



wiiA ip(x) = d Ag (x, U, V, E h ) . 

Remark 4.2. The negative exponents ofh in the upper bounds are not the optimal ones. Some care 
especially has to be taken while modifying ip or while commuting hd x with . This presentation 
is the most flexible one for our purpose. 

Proof: i) Our assumptions imply that the functions ip(x) = dAg(x> K, V, Re z) and <pu{%) = 
dA g {x, Kh, V — h, Re z), with Kh = {x £ (a, b), d(x,K)<h}, satisfy 

\(p(x) - (p h (x)\ < K atb!C h i.e. f^rj - eKa ' b ' c ' 

for some uniform constant K a ,b,c ■ Hence the function <p can be replaced by cp^ in the proof. 
Lemma lA.31 applied with a = a, f3 = b, u% = U2 = u, (f = (ph and v = e~i^u implies 

b 

2 



pb no I'D 

Re / vf> / \hv'\ 2 + I h\v 

J a J a J a 



Hence we get 



Nlk>» < J(IIMU>IHU» + ||/a|UIHU-) 



The Gagliardo-Nirenberg estimate sup xe [ Qi& ] \v(x)\ < ^-alMI^^^IMI^^)) im P lies: 

\H 2 m , h < \ (\h\\v\W\ H ^ + ^V-AmMW* 

This combined with the equivalence of ||w||^fi,h with \\he~^u'\\ L 2 + ||m||l 2 leads hnally to (|4.5[) . 
ii) We follow the ideas of [24 which consists in putting the possibly negative term of the energy 
estimate in the left hand-side. Hence the equation (Hp — z)u = f is simply rewritten 

-z)u = f + W x h u + W^u , 

and it suffices to estimate || W / 2 ' l u||_/vi [ , • The Gagliardo-Nirenberg estimate gives 

\\w 2 h u\\ Mb < c 6 _ a ||^|UJ|^|li / 2 2 IHli / 2 2 . 



25 



Applying Lemma lA.31 with a — a, j3 = b, u\ = u 2 = u, ip = leads to 















hu'\\ 2 L 2 +c||u||| 2 < 


f fu 


+ 


fw^u] 2 


+ 


[ b W 2 h \u\ 2 




J a 




J a 




J a 



Apply a second time the Gagliardo-Nirenberg estimate for 

\\u\\ m ,H < C aAc (\\f\\ L 2 + \\W?\\ L °o\\u\\ L * + hw 2 h \\ Mb \\u\\Li 



gives 



\\w 2 h u\\ Mb < ^\\w 2 h \\ Mb (\\f\\ L2 + IKIU-HU* + l\\w 2 h \\ Mb \\u\\L* 

Combined with the results of i) applied with / replaced by / + W h u, this yields: 

C 



h 1 ' 2 sup \e^u{x)\ + \\he^u'\\ L 2 + \\e^u\\ L , < ^(||/|| i2 + ||W?|k«|M| £ . 



x£[a,b] 



h 

+ IwwZWmMWl* + ^ll^ii^iiwfiu-Hi^ + ^\\w£\\ 2 Mb \\u\\ L *) 



where <ph(x) = dA g (x, if^, V, Re z) and Kh = {x € (a, 6), d(x,K' U U) < h}. With the assump- 
tions on Wi and W 2 and replacing (fh(x) by <p(x) = dA g {x, K' U U, V, Re z), we obtain (I4.6[) . □ 

4.2 Reduced boundary problem for generalized eigenfunctions 

We shall consider the boundary value problem 



(P h - z)u = f . 



hd x + <V2 

hd x - iC, 1 ' 2 



u(a) = l a 

u(b) = 4 , 



(4.7) 



with Im z and Im £ small enough w.r.t h > and specified later. Here z 1 / 2 denotes the complex 
square root with the determination argz G [ — §, if) ■ 

The case / = occurs while studying the generalized eigenfunctions of Hv,e a (0) or their variation 
w.r.t 6*o- The case l a = lb = is concerned with the resolvent estimates for the non self-adjoint 
Hamiltonians 

H% := P h , with D{H£) = |w g H 2 ((a, b)) , [hd x + i( 1/2 }u(a) = , 

[hd x - iC 1,2 ]u(b) = O} , (4.8) 
H% := P h , with D(H£) = {ue H\(a, b)) , P h u G L 2 ((a, b)) , 

[hd x + i( 1/2 }u(a) = , [hd x - i( 1/2 }u{b) = 0} . (4.9) 

Lemma 4.3. Assume V — Re z > c with ||V||£oo < i and | Im C 1 ^ 2 ! < 7~~ for Kb- a large enough 
according to b — a. Let K be a compact subset of [a, b] and set tp — dA g (x, K, V, Re z). Then any 
solution u G L 2 {{a,b)) to the boundary value problem (|4.7[) with f = satisfies: 



h 1/2 sup \e ± ^u(x)\ + \\he ± %u / \\ L 2 + \\e ±: &u\\ L 2 < 



C a ,b,c r, n i , i „ i + yl 6 ) 



a;£E [a,b] 



|4|e d 



Proof: Again the function ip is replaced by tph(x) = dA g {x, K,V — h, Re z) . Applying Lemma lA.3l 
with a = a, p = b, u\ = u 2 = u and v = e~u implies 

0> f\hv'\ 2 + f h\v\ 2 + hRc HC 1/2 ) [\v\ 2 (a) + \v\ 2 (b)] 



hRe 



.fiW 



W (a)(e ± ^4)-^(6)(e ± - 4) 



2G 



With the Gagliardo-Nirenberg estimate, we get 

\\hv'\\l 2 +h\\v\\l^2C 2 b _ a \lm C 1/2 \\\hv'\\ L A\v\\Li < C b . a \\hv'\^\\vt^h 1/2 
which implies 



141 



M4| 



(l|W||| 2 + H|| 2 ) 



1/2 Co,;,,, 



< 



hi 



141 
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provided Kb_ a is large enough according to the Gagliardo-Nirenberg constant C b - a ■ Rewriting the 
inequality with the uniform equivalence 1 1 he^ ~^u'\\ L 2 + \\e ± ~^u\\^2 with || v \\ jji.h = (\\hv'\\ 
yields the result. 

The generalized eigenfunction ?/>_ g B (k, x), fc <E K*, of Hy go (0) is the solution to 



2 

L 2 



□ 



1/2 





- fc 2 )v> = o 


inR\ - 


[a, 6} 




f ) = e~^"V(a 


") , 


?//(a + ) = e-^V'(a") 




) = e 2 


h ) , 


^(6-) = e -^y/(&+) 




■ kx 

oo,a) 


R(k)e 


. ^l( 6)+0 o)= T ( fe ) 




, = T{k)e 




^l(6,+oo)= ei *+ i? ( A; ) e " 



T{k)e l — for fc > 



This can be reformulated as the boundary value problem in (a, b) 
(P h - fc 2 )V> = in (a, 6), 

(fc > 0) 

[hd x + i(k 2 y/ 2 e- e °] i/j(a) = 2ike-^e i! f 





(k < 0) 
0, 

2ike 2 e <• 



(4.10) 



[/i9 x - i(k 2 y/ 2 e- e °] i>(b) = 

where the choice of z 1 / 2 says (fc 2 ) 1 / 2 — \k\ for k € K* . A straightforward application of Lemma l4.3l 
gives the next result. 

Proposition 4.4. Assume V — k 2 > c, || V||l°° < - and |(?o| < ch for c small enough according to 
a,b. The generalized eigenfunction ip 1 ^ o (k, ■) satisfies 



h 1 ' 2 sup le^^ffc.iJI + II^^^O'llw + lle*^^,.)!^ 

a:G[a,h] ' 

to^/i = dA g {x 1 a, V, fc 2 ) w/ien fc > , 

and with ip(x) = dA g (x,b,V, k 2 ) when k<0. 

With this first a priori estimate, the boundary value problem (|4.10[) can be rewritten 

' (P h - fc 2 )V> = in (a, 6) , 



[hd x + i{k 2 f' 2 } ^(a) = 

{ [hd x -i(k 2 )V 2 )]m 



(fc > 0) 
2%ke i> * +C( 1 ^ i ) 

d Ag (a,b,V,k 2 ) 



(fc < 0) 

dAg(">'>. t -'. fc2 ) 

g( e " l 9 o| 



Hence the difference u = i/c 9o — ip_ solves the boundary value problem 
' (P h - fc 2 )w = in (a, b) , 



[hd x + i{k 2 ) 1 / 2 ]u{a) = 

{ [hd x -i{k 2 y/ 2 )]u{b) = 

Hence Lemma 14.31 yields the next comparison result. 



(fc > 0) 

d Ag (a,b,V,k') 



0{- 



(fc < 0) 

d Aa (o.,b,V,k 2 ) 

o f ; |flo1 -), 

O(J^i). 
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Proposition 4.5. Assume V — k 2 > c, \\V\\l°° < - and \6$\ < ch for c small enough according to 
a, b. The difference of generalized eigenf unctions u = ip_ @ (k, .) — f/c (fc, .) satisfies 

,1/2 i i£M i \ i , I,, £ / £ ii , Ca. 6,c #0 

/i i/2 sup |e " u(a:)| + ||fte^ii'||i2 + ||e"«||i2 < — ^75— L , 



with 
and with 



x£ [a, 6] 

= dA g (x, a, V, k 2 ) when k>0 
<p(x) — dA g (x,b,V, k 2 ) when k<0 



Remark 4.6. With an additional regularity assumption \4@/ proves in the case 0$ = that the 
upper bound of Sup xe r Q w \e £ ^ 1 'ip h i (k,x)\ is actually 0(1) with a first order WKB approxima- 
tion. By using this result and the comparison result of Proposition \4-5\ with a bootstrap argument 
or reconsidering the complete proof of 14$ , the estimate sup x6 r a)6 i \e^~ ^ (k, x)\ — 0(1) and 

sup^guy \e !£ ^~ {4> h i q o — ip'tl )(k,x)\ — O(y^) could be obtained. Here only V € L°° is assumed 
with a possible loss in the h-exponent. 



4.3 Weighted resolvent estimates 

We complete the analysis ol the previous subsection with results concerned with the resolvent 
{Hq — z) -1 corresponding to the boundary value problem (|4.7[) with £ a = if, = . 

Proposition 4.7. Assume V — Re z > c, \\V\\l°° < ~ and | Im Q 1 ^ 2 \ < irai/i Kh- a /arye 

enough according to b — a. Let K be a compact subset of [a,b] and set (p — dAg(x, K,V,He z) . 
Then for f € L 2 ((a, b)), the function u = (H 1 ^ — z)^ 1 f satisfies 



h 1/2 sup \e ±!E ^~u(x)\ + \\he ± %u'\\ L 2 + ||e ± ^ii|| i a < 

x£ [a.b] 



C a t 



L 1 



where e ± % f = f when supp / C K . 
In particular this yields 



1 

1 P + h 



< 



Ca,b,c 



£(L 2 ) 



and z € p (jif) ■ 



Proof: Again we can replace tp by ipy l = dA g (x, K,V — /i,Re z). Lemma IA.3I with v 
implies 



\hv'\ 2 + h / M 2 + hRe HC 1 / 2 ) [\v(a) 2 \ + \v(b)\ 2 ] < \\e^ f\\ L 4v\\ } 



Absorbing the boundary term with the help of the Gagliardo-Nirenberg inequality like in the proof 
of Lemma POl and taking K&_ a large enough yields the result. □ 

Proposition 4.8. Assume V - Re z > c, \\V\\ L °° <-,\1m z 1/2 | < and I Im C 1/2 | < 
with Kb-a large enough according to b — a. Let K be a compact subset of [a, b] and set ip = 
dA g {x,K,V,Re z) . Then for f e L 2 ((a 7 b)), the difference w = (H^ — z) 



\f - (H>> - z)" 1 / 



satisfies 



h 1 ' 2 sup \e ±£ ^w(x)\ + Whe^w'U* + \\e^w\\ L , < Ca *^ ^ \\e^ 



x£[a,b] 

where e^n f = f when supp / C K . 
In particular this yields 



h 2 



/II 



L 2 , 



(H 



h z)- 1 



(Ht-z)- 1 



=1= 4- 

h 



< 



Cg, b , c \z 1/2 ~ C l/2 | 
/l 2 
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Proof: The function (H^ — z) 1 f solves (|4.7[) with £ a = it = 0. Therefore, if we set u = 
(Hi? — z)' 1 / and v = (Hq — z) -1 /, the function w = u — v verifies: 



(P h -z)w = 0. 
hd x + iz 1 / 2 
hd x - iz 1 / 2 



w (a) = -i{zV 2 - C^ 2 )v(a) . 
w(b) = i(z^ 2 - C 1/2 )v{b) . 



Then it follows from Lemma |4~31 that 



h 1 ' 2 sup \e ±S ^w{x)\ + \\he^w'\\ L , + \\e^w\\ L , < 9*&L\J-/* - f/^h 1 /* sup \e ±£ ^v{x)\ , 



x£ [a.b] 



x£ [a.b] 



and we can apply Proposition 14. 71 to the function v to get the result. 



□ 



5 Accurate analysis of resonances 

In this section, we use the approach of Hclffcr-Sjostrand relying on the introduction of a Grushin 
problem (see [30] [SB]). This section ends with a rewriting of the Fermi Golden rule (|1.2[) for the 
modified Hamitonian Hg v _ W h ■ 

5.1 Resonances 

Resonances for Hg v _ W h — Hg o v _ wh {0) are eigenvalues of Hq q v _ W h(iT~) for a suitable choice 
of t according to the resonances to be revealed. Associated eigenfunctions are the g r € L 2 (R) 
functions satisfying 

H 8 ,V-W h ( iT )9r = ZrQr , (5.1) 

with arg(z r ) G (— 2r, 0). Alternatively, f r — U-i T g r satisfies 

^9 ,V~W h f r ~ Z rfr j 

with f r e L 2 ((a + e lr IR_) U (a, b) U (b + e* T R + )) . We refer to [29] for a general comparison of 
the two approaches. Accordingly, we recover the definition of Gamow resonant functions with no 
incoming data and slowly exponentially increasing outgoing waves. 

Equivalently working with g r , the condition g r 6 L 2 (R) imposes the exponential modes in the 
exterior domain: 

ff r (x) = I g„nt(x), xe{a,b) (5.2) 



1/2. 



g~e 1 s (x a)^ ^ < a j 

where we recall that z 1 / 2 denotes the complex square root with the determination arg z£ [ — §> if )• 
According to the definition of D(Hg v _ w h (#)), this function verifies the following boundary con- 
ditions 

(ftflis - izyV r ) 3r (6+) = ^ - iz^ 2 e-°A g r (b~) = , 

and 

(hd x + iz x J 2 e iT \ g r (aT) = => (hd x + i^V 90 ) 5r (a+) = 

(with (z r e~ 26a ) = ^ 2 e~ 0a ). It follows that the interior part of the solution satisfies the non- 
linear eigenvalue problem 

H z ^ e - 2 e gr,int = z r9r,int (5-3) 

(see definition (|4.9[1 ). Conversely, given z r in the sector: arg(z r ) e (— 2r, 0) for which gy.jnt fulfilling 
(I5.3[) , it is possible to define suitable coefficients g+ and <7_ such that the function g r given by (|5.2I) 
is in D(Hg o v _ W h(9)) and solves the equation (|5.1[) . This allows to identify resonances with the 
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poles of (H^ e _ 2f)0 - z) . 

It is worthwhile to notice that this technique extends to the first Riemann sheet: in this case the 
poles of (iJ^ e _ 2 e — z) correspond to proper eigenvalues of Hq v _ W h provided that arg(z) < 4^ - 
To this concern, the following spectral characterization holds. 

Lemma 5.1. Let H. v _ W h(Q) be defined as in (|3.3[) with r £ (0, j), then 

<r P (H? T v_ wh (0))n{lm z > 0} = . 

Proof: According to the Proposition 13.61 the points in <r p [Hq v-w h (^)J ^ fi m z > 0} coincides 
with the eigenvalues of Hg v _ w h{9), in Im z > 0, for any choice of 9 with: Im 9 g (0, \ 
particular, for 9 = 9q = ir, it follows from p. lip that the operator Hjr 

this case, a p (h^ t v _ W h(iT~)j n {Im z > 0} = 0. 



In 

iT.v-w h ( iT ) is accretive - In 

□ 
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5.2 The Grushin problem for resonances 

In the previous section we got some accurate estimates for the variation w.r.t 9q of the generalized 
eigenfunctions of the filled well Hamiltonian Hq q v . Here the resonances for the full Hamiltonians 
Hq v _ W h and H^ v _ wh are considered. After reducing the problem to the interval [a, 6], we 
introduce like in 30 [20 21] the Grushin problem modelled from the Dirichlet operator with the 

z with £ = z or £ = ze 

We assume that a cluster of eigenvalues X\, . . . , A^ of the Dirichlet operator 

H h D = -h 2 A + V - W h 

exists such that 

d(\°,a(H'i)\{X h 1 ,...,X^)>c, 



potential V — W for the boundary value operator H£ 
(SU). 



according to 



c<A°< inf V(x)-c< \\V\\ L ° 

x£(a,b) 



l 

< - 

c 



(5.4) 
(5.5) 



max |A$- A°| < -h. (5.6) 

l<j<£ C 

The domain u c h will be a neighborhood of { Aj , . . . , A^ } such that 

u ch c{zeC, d(z,{Xl,...,X^})<ch} . (5.7) 

Remark 5.2. Notice that these assumptions do not forbid h-dependent X° with \X°(h) — A°| < 
since, in this case, it suffices to replace V by V — X°(h) + X° . 

Normalized eigenvectors associated with the Xj are denoted by $^ and the total spectral pro- 
jector is 



We also introduce the bounded operators 



Rn 



L\{a,b)) 



i — y Rq u 



and Rq 



L 2 ((a,6)) ->■ C 



i 

/{$>} . 
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D j =j+ q° J : D(Hp) xC ( 4 L 2 ((a, b)) x C e is invertible 

with the inverse 

H h D -z R^y 1 = (E a {z) E+ 
R+ J \ £ " Ev + {z) 

l 

E (z) = (H h D z)-\l H h ) , E+v+ = . 



Notations: We set 




£; - + (z) w + = diag(z- A£)v 



The problem — z)u = / is studied after introducing the matricial operator 



v 

where the function Xh € Co°((a, 6)) satisfies 

\\{hd x ) a Xh\\L°°{{ a ,b))<C a , a£fi and X h{x) = 1 if d(x, {a, b}) > h . 

Another cut-off function \p £ Cfi°((a,b)) will be used with a smaller support. By introducing the 
positive quantity 

S a :=d Ag ({a,b},U,V,\°), 
the cut-off ip is chosen such that for some r\ > independent of h > but to be specified later 



r o ■dd Ag (x,u,v,\°)>^ 

\ 1 i£d Ag (x,U,V,\°)< V- 



When 77 > and /i > arc small enough 

[/CC{i/iEl}c SUppV> CC {Xh = 1} • 
For z,Q £ u> c h, consider the approximate inverse 

C( } " I E^ ) ' 

where the function ph £ C(f(U~2h) satisfies 

\\(hd x ) a ph\\L°°((a,b)) <C a , a £ N and p h = 1 on Ush , 

after recalling C/ t = {.x e (a, 6), d(x, [/) < i} . In particular this implies W h (l — ph) = 0. 
A direct calculation gives 

U c {z)F c {z)=(^ B D 

with A = l+[(hd x f,p h ] {H%-z)-\l-il>)-[(hd x f,Xh]E il>, 
B = (H£ - z) Xh E+ + X hRoEo + , 
C = RtiXhE^ + (1 - Ph )(H' c l - z)-\l - $)) , 
D = RqXhEq , 
where we have used H^Xh = HjjXh = -P' l X/i • 
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Proposition 5.3. Assume the conditions (|5.4I)(|5.5D(I5.6[) and suppose z,( 6 0J c h ■ The matricial 
operator Fq{z) is an approximate inverse ofH^(z): 

U i {z)F i {z) = l + K c {z) and F c {z)H c {z) = 1 + K'^z) , (5.9) 

for h > small enough and after adjusting the parameter rj > so that \\K,q(z)\\ + \\K.'^(z)\\ < 1 
according to: 

||/C c (z)|| + \\K' c (z)\\ < C aAc e- S °' C ^' b ' cV . (5.10) 



More precisely the remainder term equals 



' {Z> I Rt{Xh-l)E^ + Rt{l-Ph){H<}-z)-\l-iP 



i> -{(hd x ) 2 , Xh ]E+ 



and is estimated by 

K c (z) 



0(e ss^) 0{^-) 

K 



,0(e 0(^) / 



Proof: Set K,q{z) = [ ^} 2 I and remember the expressions of A,B,C,D in %q{z)J : q{z) 



K21 K 2 2, 

A B^ 

C D; 

The first coefficient K\x is simply A — 1 according to the above definition. 

The coefficient K12 = B is computed by making use of H^xt — HdXIi and of the relation 
(#£ - z)E+ + R^Eq + = coming from H D (z)£ D (z) = 1 . 

The coefficient -^2,1 = C is computed after using the relation RqEo = coming from T-Ld{z)£d{z) = 
1. 

The coefficient K22 = D — 1 is computed after using RqEq = 1 . 
Estimate of K\\i For the first term, remark the identity: 

l(hd x ) 2 , Ph } = 2{hp' h ){hd x ) + {h 2 pl) , (5.11) 

where the coefficients hp' h and h 2 p'^ are uniformly bounded and supported in U211 ■ Then, owing to 
Proposition ^. 71 it is estimated with: 

|| [(hd x ) 2 ,p h ] (H{ - z)-\l - 1>)\\ < C aAc e- S °~ C ^ b,cV . (5.12) 

For the second term, we have the identity (|5.1ip . where ph is replaced by Xh and the coefficients 
hx'h and h 2 x'h are uniformly bounded and supported in {x e (a, b), d(x, {a, b}) < h} . 
By introducing a circle 70 = {V € C , \z' — A°| = -h \ , the formula 



E (z) = (H*-z)-\l-Tl h ) = -— I - -— TThs d2>. 



1 

2ut J l0 z< ~z(z>- H* ] 
and Proposition 14. 14 -ii^l imply 

So-c a , b|C „ 

\\[{hd x )\xh]E Q n < a '"' ce h3 • (5.13) 

Estimate of K21: The cut-off ( X h ~ 1) is supported in {x G (a, &), {a, 6}) < h} . Meanwhile 
we verify with the same argument as for (|5.13[) that the function u = Eoipf satisfies 

||e*«|| < CaAcim 



h 3 
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for if = d,Ag(x, supp^/i, V, A ) . The operator Rq is the finite rank operator defined by taking 
the scalar product with <F^, j = 1, ...,£. With the exponential decay of the eigenfunctions 
j = 1, ...,£, stated in Proposition 14. 11 we get: 

\\R+(1 - Xh)Eo4>\\ < C aAe e « . 

The second term is estimated like the first one of K\\ while replacing [(hd x ) 2 , pu\ with Rq(1~ ph): 

\\R+(1 - Ph )(H' c l - z)-\l - < C aAc e- S °^' b ' cV . 

Estimate of Kyi and K^i'. The operator Eq is defined by E^v + = X^-=i v j^j anc ^ the exponen- 
tial decay of the eigenfunctions j = 1, . . . , £, stated in Proposition 14. II with the relation (|5.11[) . 
where ph is replaced by Xh, yields 

\\Kvt\\ = || [(hd x ) 2 , X h]E+\\ < C a . h /—- . (5.14) 



For if 22 we use additionally the exponential decay of the $J,j' = l,...,l contained in R^ and we 
get 



_2S a 

||^ 22 || = (1 - Xh)E+\\ < C aAc ^- . (5.15) 

Left and Right inverse: When h and r\ are small enough the previous analysis says that J-"^(z)(l + 
/C(j(z)) _1 is a right-inverse of Hc.{z) and T-L^{z) is surjective. 
From the definitions (14.81) and (14.91). we have: 



H£=H h ((i^j and H% = H h ((i) . (5.16) 

After two integrations by part, we get (^H h = H h f-(C)')- With {Rq)* = R„ and with 

the notations induced by (|5.8[) and (|5.16l) . we obtain T-L((^,z) = T-L(—(()^ 7 z). The analysis 

performed to obtain (|5.10l) for /C(C^ , z) can be adapted in the case of /C(— (£)^ , z): this yields the 
surjectivity of H(— (C) 1 , z). Since 

Ker (U((* , z)) = [Ran (h{(^ , z)) *] = {0} , 

the injectivity of H((^,z) follows. □ 
Notation: When h > is small enough, we set 

The Schur complement formula 

{H* -z)- 1 = E-E + {E- + y 1 E- (5.18) 

recalls that (H* — z) is invertible if and only if the I x I square matrix E h is invertible. An 

accurate calculation of this matrix allows to identify the poles of (H^ — z)^ 1 . 

The final result comes from a higher order estimate after taking the Neumann series 

(1 + K^z))- 1 = 1 - K c (z) + JC ( (z) 2 - /C c (z) 3 + /C c (z) 4 + O(e-° S0 ^' b " V ) . 
Proposition 5.4. Assume the conditions (|5.4[)([5.5p(|5.6p and suppose z,£ € u) c h ■ Then 
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and 



E- + = - E^[(hd x ) 2 ,p h }(H^ - z)- 1 [(hd x ) 2 , X h]E+ - E^ + R+{x h - l)E+ 

_ 5Sq — C a 5 c r; 

- Ev+R+(l ~ Ph)(H% - z)- 1 [(hd x ) 2 , X h]E+ + 0(e s^) . 

(«) (") ) ' 

(2) 

= K11K12 + K12K22 : Due to the support condition when rj > is chosen small enough and 
h > is small enough, the first term equals: 

KuKu = -{(hd x ) 2 ,p h ](H£ - z)- 1 [{hd x )\ X h]E+ , 

and with the same argument as for (j5.12|) we get: 

||^n^i2|| <%^e-^, (5.19) 

where some additional exponential decay comes from the eigenfunctions appearing in Eq and the 
support of the derivatives of Xh- From the equations (|5.14[) and (|5.15[) . the second term satishes 

\\K 12 K 22 \\ < \\K 12 \\ \\K„\\ < ( -^r-> — ■ 

= K21K12 + K22 : The first term equals 

^21^12 = -B,+ (l - p h ){H* - z)- 1 [(hd x ) 2 , X h]E+ , 
and as it was done for (|5.19p . we obtain: 

II TS II ^ C a h r — 2S 

11-^21^1211 < —f^~ e " ■ 

Then, from equation (|5.15D . the second term verifies 

\\I<! 2 \\<\\K 22 \\ 2 <^fe-^. 
Estimate of K[^, , and A direct computation gives for i = 1 and 2: 

K i2 — K il K l 2 + K i2^ 22 , 

and (|5.10l) implies: 

||/c c (z)|| (||^ 2 l) || + \\k™\\) < c aAc e-^t^ + . 

Moreover, we have obtained: 

f2l , 2S 0- C a,b,c1 , t? ) , 2 s -c a b v 

K[f = 0(e R ) and = 0(e ) , 

therefore we have: 



^ 3 2 } = 0(e s^), = 0(e ^ ), 

SS 0- C a,b.c r l , . 5S 0- C a,b,c') 



Computing E h : We have 

^ + -^ + = E^il>[-K 12 +K™-K$ +K W ] + E-+[-K 22 +K™ 
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Since the operators E ip and E + are uniformly bounded, it follows from E Q ipKi 2 = and 
K 22 = 0(e-^rh- 2 ) that: 



E~+ - Eq+ = 0(e—^h- 3 ) 



and 



E~+ - E„ + = Eo^K u K l2 - E^+K 22 + E^+K 21 K 12 + 0(e~ " s °7 h ab ^ ) 

= -E-[{hd x f,p h ](H^ - z)- 1 [(hd x ) 2 , X h]E+ - E-+R+( Xh - 1)E+ 

- E^+R+(l - p h )(H* - z)- 1 [(hd x ) 2 , X h]E+ + 0( e - 5S ° . 



□ 



h 3 



5.3 Localization of the resonances 

In what follows we discuss the problem of resonances for the operator Hg g v _ w h(0)- Using (|5.18D 
and the detecting method introduced in Subsection 15.11 these coincides with the singularities of 
the matrix (i? h (z, ze~ 26 °z)) in a sector arg(z) G (— 2r, 0) for a suitable r. Here the symbol 
E \~(z, C) actually denotes the (z, £)-dependent matrix defined in (I5.17[) . 

The comparison of the Schur complements E h and Eq + stated in Proposition l5.4[ allows to state 
the following localization result on the resonances of the operator Hy_ wh 9o (0) and to estimate 
accurately their variations w.r.t 8q . 

Proposition 5.5. Assume the conditions (I5.4[)(|5.5I)(|5.6I) and fix 9 o such that \6q\ < ^ . Then for 
h > small enough, the operator H# o y_ v ^h{^) ^ as exac ^y £ resonances {zi(6q), . . . , z^(6q)} in 
LOch , possibly counted with multiplicities, with the estimate 

(_2So \ 
- 

after the proper labelling with respect to j £ {!,...,£}. 
In particular, when 

lim/i 3 e^min|A^-A^| =+oo, (5.20) 

{ 2S h" 1 

there exists T a f, c > 1, such that every disc D jth {T) = iz e C , \z - A^ l | < T^f- contains 

exactly one resonance Zj(9o) when T is fixed so that T > T a ^, c and h > is small enough. 

Proof: We look for the points where the matrix E *~(z, ze~ 26 °) is not invertible. When z £ uich, 

then \z\ < | when h is small enough and the two points z and £ = ze~ 26a = z + z(e 2e ° — 1) belong 
to uj c h- Thus, Proposition 15.41 gives 

_2S a 

\\E-+(z,ze- 28 °) - E^+(z)\\ x < M aAc ^- , (5.21) 

where the equivalent norm ||ai,||oo = max|aij| is used. 

i,j 

{_2So \ 
h e h \ 

zeC; min \z — \j \ < 2lM a ^. c — — \ and suppose z f. O/j . Then the coefficients 

I 

Er.+ of E~ + are such that for all i G {1, ...,£}, \E^ + \ > ^\ E ij + \ and E ~ + is invertible by 

i=i 

Gershgorin circle theorem. 

To conclude the proof, we have to compare the number of resonances to the number of Dirichlet 
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eigenvalues in each connected component £lj,h of {Qj,h — &j',h is not forbidden). Defining 
E-+(t) = Eq + +t(E~ + -Eq + ) for < t < 1, the number iV(t) of points in Q jih such that E~ + (t) 
is not invertible, is constant on [0, 1]. Actually, note first that (|5.21l) implies that for all f £ [0, 1], 
E '"(i) is invertible when z £ dflj^hi using an argument similar to the one used for E h outside 
ilh- Therefore, for any to £ [0, 1] the analyticity of E '"(to) with respect to z implies 

inf | det £7" + (t )| > , 

and for S small enough, the estimate: 

|dct£T + (t + (5)-dct£r + (t )| = \dct(E-+(t ) + S(E- + -E Q : + ))-dctE-+(t Q )\ = \5\\R{t ,6)\ 

< inf |dct£T + (t )| < | det J5 _+ (t Q )| 

z£d£l jih 

holds for all z £ dflj^h after noticing that the function |i?(to,<5)| is a bounded polynomial of to, 
S and of the coefficients of Eq + and E h . The functions det(E '"(to)) and detE '"(to + 8) are 
holomorphic functions of z £ ujch such that sup z€ Q n . h I d dcfg-+fo)^ — 1| < !• Thus, Rouche's 
theorem implies N(to + 8) = N(to). The function N(t) is continuous on [0, 1] with integer values. 
It is constant. 

Assuming e~ h 3 \Xj — Xj,\ — > +oo for all pair of distinct j,f, implies Clj^ C Djj l (R) for all the 
j's with Dj : h(R) n Dji t h(R) = if j ^ f when R > 2iM a ^ c and h is small enough. This yields 
the last statement. □ 

Remark 5.6. In the above proposition the term resonances is used for the eigenvalues of the 
operator H _ 2 e > which in principle may still have a positive imaginary part. In the particular 
case of 8q = it. t £ (0, ?), the result of Lemma \5.1\ implies that these eigenvalues must lay in 
the lower half complex plane. On the other hand, the result of next proposition and the lower 
bound on |lm Zj(0)| (see Proposition 1 5. 8}) implicitely yields: Im 2?(#o) < on a suitable range 
of \8q\. Under each of such conditions the points Zj(8o) corresponds to resonances of the operator 
Hg o v _ W h(0) as defined in Provosition \3. b\ 

The next Proposition localizes the resonances z^(9q) of Hg o v _ W h(Q) with respect to the reso- 
nances z'j := Zj(0) of Hq V _ w h (0) by making use of the comparison between E~ + (z, ze~ 29 °) and 
E-+(z,z). 

Proposition 5.7. Assume the conditions fO j) (|B1 ]) ([5^6 ) and < |0 O | < ^ ■ Then for h > 

small enough, the matrices E h of Proposition \5.4\ associated with ( = z and £ = ze~ 26 satisfy 

(_2Sq \ 
/ 

If additionally (|5.20j) is assumed the variation of the resonances around X° fore-T% < \9 \ < M 
and 9q — is estimated by 

max |«?((9o)-*il 
je{i,...,e} 

Proof: For z £ ujah and 9q such that \9q\ < the Proposition 15.41 implies that: 

E-+(z, ze- 29 °) - E~+(z, z) =Eo[(hd x ) 2 ,p h ]D[(hd x ) 2 ,Xh]E£ 

b.cl. 




+ E^+R+(l - Ph )D[(hd x ) 2 ,Xh]E+ + 0(e 

5S 0~ C a,b 

=I + II + 0{e n ), 
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where D = (H^ — z) 1 — (H^ e _ 20o — z) 1 . The operator E Q being bounded, the first term / is 
estimated as we did for (15 . 19[) where (H^ — z)" 1 is replaced by D and we use Proposition 14.8 
instead of Proposition 14.71 This leads to: 



2 Sp 

e h 



2S 



For //, using the exponential decay given by the operator Rq , we get: 



|//|| <C a , 6 , c |0 o | 



h 3 



The assumption e~ < 1 6*o | ensures that the remainder O (e 27r^ ) is absorbed by | 6q |e s 11 /i~ 3 
as h ->■ . We have proved (|5.22l) . 

When (I5.20[) is verified, Proposition 15.51 savs that every disc Dj t h(T) — {z e C, |z — A^ l | < 
Te-^/i- 3 } for 

any > T a jj C , contains exactly one resonance z^{9q), and in particular one 
resonance Zj when 6q = 0. Hence the matrix E ^(z,z) has only simple poles and its inverse is 
the meromorphic function 



The matrix Aj is nothing but the residue 

while the function F h (z) is a holomorphic function estimated via the maximum principle by 



sup \F h {z)\ < sup 

4 

The estimate of Proposition [ITU says 



~4~ 



|( J B-+(z,z))- 1 |+^ 



14 



i i z ~ z i 



(5.24) 



|^-+(z,z)-S -+(z)| < C aA 



e h 



while we know \(E + (z)) 1 | < maxi<j<£ \z — 1 . After writing 

{E-+)- 1 = [1 + (E^ + )- 1 (E-+ E "+)] 1 (S -+)- 1 
we get for T > max{T aif , !C , 2C aj6 , c } 



(5.25) 



sup z))- 1 ! < 

zedD j:h (T) 

and finally the uniform bound for the residues 



< 



t[i - 



T 



max lA^I < 2. 

i<j<e J 

The holomorphic part F h (z) is then estimated with (I5.24p . Actually, the first term is estimated 
with the help of (|5.25l) while the second term is treated with the above estimate of A 1 !- and by 



making use of maxi^x^ \zj - Aj| < Ta^^ 



sup \F h (z)\ <% 
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For all z € ujch \ {zj 1 , . . . , z^} , the inverse of E + (z, z) is thus estimated by 

2 C'_ hr 



\{E-+{z,z))- l \<Y, 



3=1 .*"*j fc l 



We now write for z ^ Wd, \ {zj 1 , . . . , z^ } 

£T+(z, ze- 2e °) = E~+(z, z) [1 + {E-+(z, z))" 1 (£T+ (z, ze" 200 ) - £T + (*, z))] 
Due to the estimate (|5.22[) the condition 



3\ 



mm z — z„- > i - 



je{i,...,£} 3 h 3 
implies 



\(E-+(z, z))-\E-+(z, ze- 29 °) - E~+(z, z))\ < 



h 3 ^ 

+ c' aAc h- x 



T\e \ 



where the right-hand side is smaller than 1 if T > 4£C" b and ft. > is small enough. Outside 



! |z e C, |z — z^\ < T\6o\e ~~h~h 3 \,E + (z,ze 2e " ) is invertible. For such a T we have proved 



max |z' l ((9 ) - z' l | < T 



e"— V 



J J h 3 

□ 

5.4 A Fermi-Golden rule 

In [5T] , a Fermi Golden rule for the imaginary parts of resonances Tj = — Im zj l in the case 
8q = has been introduced. It plays a major role in the analysis of the nonlinear effects studied in 
[15] [20 21 [45] [2|)] [IB] f° r it expresses accurately how the tunnel effect between the resonant state 
and the incoming waves is balanced between the left and right-hand sides. By assuming 

lime^ min |A^-A?,| = +oo with Su < , (5.26) 

h^o i<j<j'<e 8 

which is stronger than (|5.20[) , the energy range of A € K associated with the resonance z^ is 

given by \ \ — Zj\ < e~~^ . When ^ (±\/A, .) denote the generalized eigenfunctions of the filled 
well Hamiltonian Hq V (0) at energy A defined in section l4~2l and j £ {1, ...,£}, denote the 
normalized eigenfunctions of the Dirichlet Hamiltonian given in f|4.3[) . the formula 

r, + o(r, ) - ^-j= > — , (5.27) 

for all A G R such that |A — Zj\ < e~~^ , has been proved under additional assumptions about the 
localization of the $^ within supp W h . We refer to Proposition 7.9 in [21] and to the subsequent 
explicit computations in Sections 7 and Section 8 of |21) for the details and in particular for the 
lower bound. 

We shall assume that the formula (|5.27[) is true when 9q = and check that it remains true 
when 6>o ^ is small enough. We shall use the notation 



I?£ = - Im z) and T^(9 ) = - Im z%(9 ) . 
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Proposition 5.8. Assume the conditions (|5.4p(l5.5[)(|5.6p(j5T2l)))(|5.27p and take 9o such that e < 
|$o I < then the Fermi Golden Rule 



I^(flb) = — r ^ nVV '" V, ^ / ' 4 ^' ^ — ^' +0 (T h 1 )+0[\e G \ e - J f \ (5.28) 



ZioZcfe /or aZZ A € M smc/i £/ia£ | A — z' 1 1 < e ^ 



j 

In particular, when limft_j.o h~ 5 Oo — we have 



2S 

e 



^ )+ ^ W ). |(w ^-*<^-'^ >i ' +i <r* > -* ( " A -'-* ?>1 ' ^ 

Proof: Proposition 15.71 gives: 

_ 2So 

|r£(0 o ) - rj| < c Oi6iC |0 o |L_ . (5.30) 

Let /(6» ) = (W h tp^ 6o (VX, .), The pointwise and L 2 weighted estimates of u = ^^{VX, .) - 
Vc o(V% •) stated in Proposition 14.51 with tp(x) = dA g (x,a, V, A) say 

,1/2 , !£M. i m , n £ I, , Ca,b,cl^o| 

ft ' sup |e h + ||e h u||L2 < 



xe[a,6] h 3 / 2 



while Proposition 14.41 gives 



ft 1 / 2 sup | e -^^ 9 (VA,x)| + ||e^ )flo (VA)|U 3 <-^, 



with the same estimate for 6$ = . Moreover, the exponential decay of stated in Proposition ^. II 
can be written as 

ft 1 / 2 sup |e^^(x)| + ||e^$^|| L 2 < , 

xe[a,b] ft 

with = d J 4 g (x, {/, V, A^). Recalling that 

\\W^\\ L ^<-, \\W£\\ Mb <~, suppW h C {d(x,U) < ft} , 
c c 

and So = t*A fl (tf, {a, 6} , V, A ) < d Ag (x, {a, b} ,V,X) + 0{h) when x e U h and |A - z h } \ < . 
Hence we get with our assumptions 

\f(0o)-f(o)\<c" b je x 



h 3/2 

We obtain 



and \f{6 )\<C'l b 
\f(0o)\ 2 |/(0)| 



h 5/2 

e b 



4ft\/A 4/iVA 



ft 5 



□ 



Remark 5.9. In /£!/, i/ie Fermi Golden Rule (|5.27[) ftas been studied with W h € L°°((a, &)). 
Nevertheless it can be proved in cases when the singular part does not vanish by a direct 
analysis like for example when W^ 1 = hS c . The presentation of Proposition \5.8\ shows that the 
stability result w.r.t to 6$ holds in this more general framework and leaves the possibility of further 
applications. 



39 



6 Accurate resolvent estimates for the whole space problem 

In the previous sections E] and [SJ we got accurate resolvent estimates with respect to h > for the 
problem reduced to the interval (a, b). We use here this information in order to derive accurate 
resolvent estimates for {Hg o V _ W h ($o) — z)^ 1 when 9 — ih N " , N > 1, which are essential in the 
justification of the adiabatic evolution. 

6.1 Localization of the spectrum 

The results of Corollary 13.41 Proposition 13. 61 Proposition 15.51 and section |5"7T1 can be summarized 
with the corresponding assumptions. 

Proposition 6.1. Assume that V € i°°((a,6);M) and W h = W{ 1 + W% € M b ((a,b)) are real 
valued with the hypothesis (|4.1j) . 

d {a , b) <V , \\V\\ L «<- , \\Wt\\<- , ||W^||<-, 

c c c 

Wi and W£ are supported in the domain Uh = {i£ (a, b); d{x, U) < h} where U is a fixed compact 
subset of (a, b) . Assume also do — ih N ° with No > 1 and h < ho . Then: 

a) cr ess (H do V _ wh (6o)) = cr ess (H ND V _ wh (6 )) = e" 2 °M. + ; 

b) The equality p.34j) written with 9 = 9o 

( H 9o,v-w h ( n) - z ) = ( H ND,v-w h ( e o) - z ) 

4 

- Yl (Bq(z,9 ,V -W h )- A)' 1 B jk (j(e k ,z,9~o),-) mm l(e l ,z,9 ) (6.1) 

holds as an identity of meromorphic functions on C \ e~ 28l) K + . 

c) In S/j = {z6C,c< Re z < inf X £(a,b) V(x) — c, | argz| < 2| Im 9q\} , the (discrete) spectrum of 

Hq o v _ w h(&o) is made of eigenvalues z h which satisfy Ker(_ff^ he _ 2 „ — z h ) ^ {0} where 
is the operator defined in (14. 9j) . When the spectrum of the Dirichlet Hamiltonian tr(iJ^) is 
made of clusters such that (|5.4p and (|5.6p are satisfied, then there exists K a ,b,c > such that 

Vz h € E h n <r(H^ v _ wh (e )) , d(z\ a{H h D )) < . 

The resolvent of Hg v _ W h{&o) will be studied in a domain surrounding a single cluster of 
eigenvalues z h , i.e. with lim^o * = with a distance to o-(H' e l v _ W h(fio)) bounded from below 

Definition 6.2. The complex domain Gh(X°) is chosen accordingly to the constants a, 6, c involved 
in the assumptions on V, W h , the hypotheses (|4.ip(|5.4l)(|5.5p(|5.6p and to 9 = ih N ° , N > 1: 

( h N ° 1 
G h (X°):= IzeC, | Re z - A°| < E aAc h, \aigz\<h N ° and d(z,a(H^))>- \, (6.2) 

for some constant S a! b jC > chosen large enough. 
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6.2 Estimates of the finite rank part 

We study here the finite rank part of (|6.1[) : 

4 

T h (z,V-W h )= J2 {Bq(z,9 Q ,V -W h ) - A)^ 1 B jk (j(e k ,z,e ),-) L2ml (e l ,z,e ) , (6.3) 

i,i,fe=l 

and its variations between the case W h — and W h ^ . Every factor will be considered 
separately. Hence it is convenient to keep the notation jv-w 1 ' {§.%, z i 0) f° r the total potential 
V - W h and 7y(e i; z, 9) when W h = . 

Proposition 6.3. Assume the hypotheses P~Tj) ((5^ (153)) (pTBl) . Take 9 = ih N ° , N > 1, and let 

G h (X°) be the set defined by ([S2|). The matrices (Bq(z,9 a ,V-W h )- A)- 1 and (Bq(z,9 a ,V)-A)- 1 
verify the uniform estimate 

Vz e G h (X°) , \(Bq(z, 9 , V — W h ) - A)~ 1 \ + \(Bq(z, 9 , V) - < C a>b<c h , (6.4) 

in any fixed matricial norm. 

The functions j\> satisfy the uniform estimates 

supp7v(e 2:3 ,z,#o) = supp7v(e 2i3 ,2,6io) C [a, b], (6.5) 
supp7v(e 1 ,z,6» ) = supp7 V (e l7 z, 9 ) C [6,+oo), (6.6) 
supp7v(e 4 ,2:,6»o) = supp7v(e 4 ,z, 9 ) C (-co, a] , (6.7) 

max||7 V (e l ,z,6»o)|| H i,h ((a ^ )) + Htv^,*, 0o)||ffi.h(( o ,&)) < "pyr ' ( 6 - 8 ) 

max||7v(s i ,2;,0o)||iri,f(R\[ o ,6]) + Wlvi^, z ,9 a )\\ H i, h(maM) < h(N ^ )/2 (6-9) 
/joZds when z € G h (A°) and V = 7 - VF' 1 or V = V . 

Moreover the differences (Bq(z, 9o, V — W ) — A)^ 1 — (Bq(z, 9q, V) — A) -1 and r ) v _ W h — -fy are- 
estimated by 

\(Bq(z, 9 0) V-W h )- A)- 1 - (Bq(z,9 0) V) - A)~ 1 \ < C aAc e~™ , (6.10) 

^flhv-w^{§.i,z,9 ) - 7v(e l ,2,6 l o)|| H iA(( a: f,)) < G a , 6jC e s* , (6.11) 

with S = d Ag ({a,b},U,V,X°), 

lv-wh{e,i,z,9 ) = lv(e t ,z,9 Q ) , /or i = l,4, (6.12) 
when z € G/j(A°), with the same result for (z,9q) . 

The matrix Bq(z, 9q, V) — A is expressed with boundary values and we will need the next lemma. 
Lemma 6.4. Assume (|5.5[) . 9o = ih Na , No > 1, and set 

ip b (x) = d Ag (x, b, V, A ) , (p a (x) = d Ag (x,a, V, A ) , 5 = ^({a, 6} , U, V, A ) . 
For any z € G/j , £/ie solution u 2 ( resp U3 ) to 

(P h - z)t2 2i3 = {-h 2 A + V — z)u 2t3 = 
M 2 (a) = 0, u 2 (o) = 1 (resp. u 3 (a) = 1, u 3 (6) = 0) , 

verifies: 

Zk , -/ 



resp. 



M|ffi,h((a,b)) < C a ,b,ch 1/2 , || e ^M2||l,2(( a ,6)) + ||e ^M 2 | 



1^(6)1 <%^, irm^l^Ca,^- 1 , |«i(o)|< C '" — 



/l 2 



resp. \u' 3 (a)\ < , | Im u 3 (a)| < a,,^" " 1 , |t&(6)| < 



^ G a ,b.c 
~ h 1 / 2 


, (6.13) 


^ C a .b,c 


,(6.14) 


2 S 


(6.15) 


h 


(6.16) 
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For any z £ G/ l (A°), the solution U2 (resp. 113) to (|3.24[) rewritten as 

(P h - z)u 2 ,3 = (-h 2 A + V - W h - z)u 2 3 = 
u 2 (a) = 0, u 2 (b) = l (resp. u 3 (a) = 1, u 3 (b) = 0) , 

can be compared with 1L2 ( resp. u 3 ) according to 

max \\ui - Ui\\ H i,h + \u'Aa) - u'Aa)\ + {u'Ab) - < C a ^_ c e~~^ 

*e{2,3} 



(6.17) 



Proof: Let us first focus on the case W h = 0. It suffices to study the case of U2, the result for u 3 
being deduced by symmetry 

Consider a real valued function uq g C°° such that uo(b) = 1, suppuo H [a, b] C [b — /i, 6] and: 

\\(hd x ) a u \\<C a , a EN, 

and set u 2 = u n + v where v solves (Hp — z)v = f with / = (hd x ) 2 u n — (V — z)u$ . We keep the 
notation (|4.2p for the Dirichlet Hamiltonian associated with P h . 

Owing to V — Re z > c, the variational formulation of (Hp — z)u = / with || /\\l 2 < Ch 1 ^ 2 provides 
< Ch 1 ! 2 which yields the first estimate of (|6.13l) . With the equation —h 2 v" = — (V—z)v+f 
and applying Lemma [A. II to v(b — x), we get < ^%jt an d hence the first estimate of (16. 15[) . 

The second estimate of (|6. 15|) comes similarly from | Im v'(b)\ < Ch N °~ 1 . It suffices to write that 
w = v — v solves 

(H D — z)w = (z — z)uq + (z — z)v = g , 

where the right-hand side is estimated by \\g\\i, 2 < C/z 1 / 2 1 Im z\ < Ch No+1 / 2 . The estimate for 

follows the same arguments as for |5'(6)| with h Na+1 / 2 instead of h 1 ^ 2 . 
The second estimate of (16.131) is a direct application of Proposition 14. ll -i) to (H D — z)v — f while 
noticing that d,A g (x, b, V, A ) can take the place of dA g (x, supp /, V, Re z) because | Re z — \°\ = 
0(h), supp / C suppuo C [b — h, b] and d.Ag(-, V, A ) is uniformly Lipschitz on [a, b] 2 . On the in- 
terval [a,a+h] the second derivative of v satisfies —h 2 v" — — (V— z)v so that ||{t2(a+-)||.H" 2 .'>((o,/i)) < 

Ce 



jjj2 — and we apply again Lemma IA.1I 

The difference v = 112 — U2 solves the Dirichlet problem 

(P h - z)v = W h u 2 
v(a) = v(b) =0, 

which means v = (H D — z)~ 1 (W h v,2)- It suffices to apply Lemma lA.21 with || W"'U2||#--i < 

r" -4? 

C a , b \\W h u 2 \\M b < C aAc e-% . With z € G h (X°), this gives \\v\\ h i,h < %^r- ■ In [a.a + h] 
or in [b — h,b], the equation for v is simply h 2 v" = (V — z)v and we use again Lemma [A. II to 
conclude . □ 
Proof of Proposition l6T3l 

a) First consider the estimate of (Bq — A)^ 1 in (|6.4|) and its variation in f|6 . 10|) . The explicit 
form of the matrix (Bq(z, 6 , V - W h ) - A), using the definitions of A,B and q(z,9o,V) given 
respectively in p.33p and (|3.26p . can be written 



with 



M(z) = - 



_36 a 



ihe" a 
Vze 2e 



'2 

fa 

e 2 



(Bq(z,6 ,V-W h )-A) 
\ 



M(z) + R(z) . 



u' 2 (b) 



-e 2 



u' 3 (a) 



ihe"o 

_3flo 

— e 2 



R( Z ) - Y 2 



u' 3 (b) 0\ 










-u' 2 (a) 
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where w 2 ,3 has to be replaced by u 2 ,3 when W = . According to f|6. 15[> f|6. 16[> and (|6.17p 



\\R{z)\\ < C aAc e™ . 
The inverse matrix (Bq(z, 9q, V — W h ) — A) is formally given by: 

(Bq{z, 9 Q , V-W h )- A)' 1 = M(z)- 1 (1 + R(z)M( z y 1 y 1 



An explicit computation gives 



M(z)- 1 = h 2 



1 

A+ 



iSL 

e 2 



u' 2 {b) 



Vze 2e o 



1 



_3£q 

-e 2 
-ti' 3 (a) 



ihe"o 

— e 2 



with 



A- 
A" 



«' 2 (i>), 



V ze 2e o 
9 ihe e o „,/ 



(6.18) 



(6.19) 



Using |Im u 2 (b)| + | Im u' 3 (a)\ < C atbtC h N °~ l due to (|BTi^ ([6"T5j) and ([637)1 . leads to the lower 
bounds for A + and A - . With our choice of the branch cut: J 1 ° 8q = ±^75, depending on argz, 



and one has 



|A+| > Re A+ > cos(Im O ) - r-r |lm u 2 ( b ) Re v^-Re u' 2 (6)Im y/z\ 

\ z \ 

>cos(Im «o)-C|^ (h N °- l y/\z\ + ^\Tm V^l 

For z G Gh, we have | Im < Ch N ° and one finally gets: |A + | > i. The same lower bound 
holds for A". Thus M(z) is invertible with: M(z)^ 1 = O (h) in any fixed matrix norm. From 
relations (|5T8"f and ([rJT9l . it follows: 



(B<z(z,fl ,^-VF' l )-A) 1 = M(z)- 1 I 1 + 



' 3S o 

e is 



(6.20) 



which is a rewriting of (|6.4I) . The estimate (I6.10|) of the difference is due to the exponentially small 
size of \u' 2 ^(a, b) — u' 2 3 (a, b)\ stated in (|6.17l) . 

b) We shall now consider the estimates for 7v(ej) with i = 2,3, V = V or V = V — W h . Actually 
it suffices to remember the equations (13.231) , (I3.24[) and the definition of the coefficients (|3.25p 



lv-w h (s.i, z, 6> ) = CiUi — ±^2 u i ' i = % 3, 

where the functions 1*2,3 do not depend on 9o and are given by (|3.24p . It suffices to use (16. 13)) for 
(|6"T5j) and (|6~T7|) for (jfjTTI) . Changing O into ^ has no effect and z G Gh(A°) when z G G h (A°) . 
c) The functions 7v(£i 4j £7 #o) do not depend on the potential V: 



7tei,2,0o) 



!(6,+oo)e h , 7(e 4 ,^,6'o) = — ^=l(_oo,a) e " 



20, , 



hy/ze 200 



from which (|6.12[) follows while (|6.9D comes from 

C 



h(e.i,z,e )\\ H1 , h{ma 



b]) 



< 



him (VI^) ~ h N < 



-Il|+T' when z G G/,.(A°) . 



□ 
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6.3 Resolvent estimates 

We gather the information given by the Krein formula (|6.1[) and the control of the finite rank part 
T h (z, V) given by Proposition! 1 



Proposition 6.5. Assume the hypotheses P~Tj) ([53 ]) ([S"lj ]) . Take 6 = ih N <> , N > 1, and let 
G h (X°) be the set defined by (|6~2j) . 

a) For V — V or V = V — W h , the resolvent (-ffe .y(#o) — z)^ 1 is estimated by 

Vz g G/j(A°) , ||«,v(0 O ) - 2)- 1 |U(£»(R) ; ffi(R\{.,6})) < ^| , (6-21) 

G & 

Vz G G/^A ) , ||(-ff e0i y_ H /''( 6 'o) ~ ^) _1 Xll£(_H--i((a,b));ifi(R\{i,6}) ^ hNo+3 ' ( 6 - 22 ) 

/or any /ixed X € Cg°((a,6)) . 

b) The difference of the resolvents equals 

(Hg Q V _ wh (9 ) - zy 1 - (Hg o V (9 ) - zy 1 = (H ND V _ wh (6o) - zy 1 

-(K Dy (9o)-z)- 1 +RT h (z) 7 (6.23) 

with 

Vz G G h (X°) , \\RT h (z)\\ c(L 2 {m . m{maM)} < C aAc e~^ , (6.24) 
Vz G G^A ) , ||i?T' i (z)x|| £ (H-i(( a , 6 )) ;ffl ( R \ {a , 6} )) < C aAc e~^ , (6.25) 
/or any fixed \ € Cg°((a, 6)) . 

Proof: a) The formula flU]) says for V = V or V = V - W h : 

(KA^-zy 1 = (H h NDy {e )- z y l 



4 

- J2 (Bq(z,6 ,V)- A)". 1 B jk ( 1 {e k ,z,e ),-) L2m 7(^,2,(90) 

= « D ,v(^)-^) _1 -T' l (z,V). 
By Proposition 1631 actually by ([0 )) (|6l? l) and (gTS}, the term T h (z, V) satisfies with 7V > 1 

II I l«, V)||£(L2(R); ff i(R\ {a ,h})) < . TTTv-TT ,Ar n/2+2 ,31 X T < 



MTV-* VWII <T Ca,b,c,xh 1 s Ca,b,c,x 

l|T («,V)xlU(ff-i((«,6))iffl(H\{a,6})) < + X fc2 - "^^+3~ ' 

for any fixed x € Cg°((a, 6)) . 

It remains to estimate the first term. The worst case is for V = V — W h since G/ t (A°) lies 
around elements of a(H% = -h 2 A D + V - W h ): 

According to Lemma E21 with ||/||if-i."((a,i>)) < jMIWh- 1 and ||tt||.ffi((a,5)) < \ IM|hi>(( ,6))) the 
inequality 



£(ff- 1 ((a,6)),.H3((a,5))) ^ 2 \d{z, o-(ITn)) / /i Ar °+ 2 
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holds for z e G h (X a ) . 

The resolvent of the Neumann Laplacian ^— h 2 A^ a ^ — (j can be written 

(-^^R\ t „6] - C) 1 = (-h 2 A^ [a>b] + l) 1 [l + (1 + C) (-h 2 A^ [aM - C) _1 



with 



2 A N 



-h z A 



R\[a,b] 



|£(L 2 (B\[a,b]);Hi(B\[a,6])) 



< 



it is estimated by 



II (-^A^ M] -C 



I £(L 2 (R\[a, b^-H 1 (K\[o,6])) 



< 



1 + 1±K 

d((,R + ) 



For 2 e G h {\°) and C = ze 29 °, the distance d(ze 
we get 



20,. 



II (-^^[0,6] -^°°J ||£(L 2 (R\[a,b]);Hi(R\[a,6])) < ^JVo+T " 

Putting all together gives (|OT|) and (j6T22l) . 

b) For the difference of resolvents, it suffices to notice that 

RT h {z) = T h {z, V) - T h {z, V - W h ) . 

Hence it is the difference of trilinear quantities of which every factor is estimated by ji l No with 

c h s ° 

variations bounded by a ' h e m, . This ends the proof. □ 



„20 r 



), is bounded from below by Ch N °. Hence 



Co,, 



7 Adiabatic evolution 

We consider now a time-dependent real valued potential V{t) — W h (t) = V(t) — W^(t) — W^it) 
supported in [a, b] with 

Ah M 2 

W?{x,t) = w n C~ X h n '\ t), W£(t) = (* h (t)h6{x-x j2t2 ), a h (t) > , (7.1) 



i2=i 



where the Xj's are fixed (independent of h) distinct points of (a, b) and the supports suppiu^ are 
contained in a fixed compact set. The functions V(.,t), ti;^ (.,£), otj 3 (t) are possibly /i-dependent 
C K functions with a uniform control of the derivatives and which impose a uniform control of 
(I4.1[)(|5.4[)(|5.5[)(|5.6[) . Namely we assume that for some X°(t) the estimates 

t max r] \\%V(t)\\ L oo + ||0 t %i(*)IU- + \d k t a, 2 {t)\ + \d?X°(t)\ < ~ , (7.2) 

< k < K 
< 31 < M- L 
< 31 < M 2 



Vx€[a,b], c<X°(t)<V(x,t)-c, 



(7.3) 



hold for all t € [0, T] . Actually, the regularity of X°(t) can be deduced from the other assumptions 
possibly by replacing it initial guess by the mean energy value iTr [ifg (i)IIo(t)] with Ilo(i) = 

This X°(t) is moreover assumed to be the center of a cluster of eigenvalues of the Dirichlet 
Hamiltonian H^(t) = -h 2 A D + V(t) - W h {t) on (a, b): There exist A?(t), . . . , X^(t) € a(H%(t)) 
such that 



d(X°(t),a(HUt))\{*i(t),...,X h e (t)})>c, 

max |A^(t) — A°(t)l < - . 
i<j'<^ J c 



(7.4) 
(7.5) 
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The operator Hq q v(t)-w h (t)(®°) ^ s studied here with 

8 Q = ih N °, N >1. 

According to Proposition IHH] and Definition ^. 21 the complex domain Gh(Xo(t)) surrounds £ eigen- 
values z\ (t), . . . , Zg(t) of Hg g v(t)-w h (t)(® } anc ^ ^ s distance to the spectrum remains uniformly 
bounded from below 

■h N °. 



min d{G h (X (t)),a(H e(uV(t) _ wh{t) (9 ))) > — 

te[0,T\ <^a,b,c 

The spectral projection associated with the cluster of eigenvalues {zj 1 , . . . , z^ 1 } is given by 



2in 



r h (t) 



(z - B(9 ,v(t)-w'>(t)(^o)) 1 dz , 



(7.6) 



where T h (t) is a contour contained in Gh(X°(t)) . 

When K > 1, the parallel transport $o(t, s), t, s € [0, T], associated with (i"b(i))te[o,T] is given by 



d t $o + [Po, d t P ] $ = 



(7.7) 



is well defined and satisfies 

VM€[0,T], flj(t)$o(*,») = *o(*,*)fl)(»). 
The time-scale is given by the parameter 

e = e~~ , with r > fixed. 
When the assumed regularity is large enough, K > 2, Proposition 13 . 71 -d ) ■ the Cauchy problem 



ied t u = H$ oV(t) _ wh(t) (6 )u, t>s. 

u(t = s) = 1A S 



defines a dynamical system II s (t, s), < s < t < T, of contractions on 1? 



(7.8) 



Theorem 7.1. Assume ([712 ]) ([Til ]) ([7T3" ]) ([775 ]) wt/i if > 2 and tafce # = i^ , No > 1 , £ = e~£, 
t > . Lei Po(i) &e i/ie spectral projection (|7.6p . Zet r belong to C°([0, T]; L 2 (R)) and fei r s belong 
to L 2 (R). for s G [0, T] tafce an initial data u s € L 2 (K) smc/i i/iai Po(s)u s = u s . 
Then the solutions u h and v h to the Cauchy problems 



ied t u h = H$ 
u h {t = s) = u s + r 



■"' = ■ a e ,v{t)-W(t)( 9 o) u + K*)> *> s > 



(7.9) 



satisfy 



ied t v h = ^ (s,t)Po(t)(Hl !nt) _ wh{t) (e Q ))Pomo(t>s)v h , t>s 
v h (t = s) = u s , 



max \\u h (t)-$ (t,s)v h (t)\\ < C aAc . T . T .s 

t£[s,T] 



e'-'Wus 



— max ||r(t)|| 
e te[s,T\ 



(7.10) 



Remark 7.2. The estimate with the source term r(t) can be improved if P(t)r(t) = r(t) after 
reconsidering the proof of Corollary \B.2\ in the Appendix (possibly with a higher order starting 
approximation with K < 2). Nevertheless the accuracy of the result may depend on the assumptions 
for r(t). We prefer to postpone this kind of improvement to a subsequent work when hypotheses 
for the source term are naturally introduced. 
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Proof: a) When UQ Q (t) denotes the solution to (|7.8p associated with r s = and r = 0, the 
contraction property of U e (t, s) implies 

max \\u{t) - u 00 (t)\\ < \\r„\\ + \ max \\r(t)\\ . 
te[s,r] e te[s,T] 

Hence, we can forget the remainder terms and simply prove the estimate (|7.10[) when r s = and 
r = 0. 

b) We consider the operator A e (t) — ^ (He ,v(t)-w h (t)(@o) — X°(t)) and we notice that the domain 
j(G h (X°(t)) - A°(t)) = | \{z- A°(t)) , z e Gfc(A (t))} 

contains the contour 

r e = 1 {r h (t) - A°(t)} , 

which can be chosen independent of t € [0,T]. Then the projection Po(t) is nothing but 

Po(t) = ^~ I (z-A^t))- 1 dt. 

Hence it suffices to verify the estimates of dj? (z — for k < K + 1 and £ € [0, T] in order to 

apply Theorem IB . 1 1 and additionally the uniform boundedness of ||Po(*)|| an d ll^t-fo(*)|| hi order 
to use its Corollary IB. 21 

Like in Appendix [Bj we use the notation g(e) = 0{e ) in order to summarize 

V5>0,3C g , s >0, \g(e)\<C g , s e N - 5 . 
For z £ T e , the fc-th derivative of [z — A £ (t))~ 1 has the form 

- A^t))- 1 = ^--Jm (* ~ A*)-^ {-iV + l W h )}(z - A 6 )^ 1 . . . 

Jl + ■ ■ ■ 3m = h 

U > 1 

...[d( m (-iV + iW h )](z- A*)' 1 , 

where the numbers Cj lt ...j e are universal coefficients. 
Remember 

the support condition 



||^V(t)|U (ia(R)) <\. 



cP t W(t) = X (x) [dlW{t)\ x(x) 

with x € C^°((a, £>)), which entails 

\\d J tW{t)\\ c(H i i{aib)) . H -i {iam < -. 
Hence the resolvent estimates of Proposition 16.51 imply 

max \\d?(z - A^t))- 1 ]] < - ^f'" = <5( £ °) • (7.11) 



Meanwhile the length |r s | is bounded by 0(1) and therefore the conclusions of Theorem IB. II are 
valid. 

Now comes the final points, which are the uniform boundedness of ||Po(*)|| an d ||<9tPo(£)ll; m order 
to refer to the more accurate version of Corollary IB. 21 
c) For Pq (t) , we write 

Po(t) = ^J r ( z ~ ^(t))- 1 dz = ^- J (z - ffflo.vW-v^wW)" 1 dz , 
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and we use the formula (|6.23p in the form 

(K,v-w»(0o) - zT 1 (H h D z)- 1 = e M0 (-/i a A^ [O)6] - ze 26 ")' 1 

+ «,v(0o) ~ z)- 1 - {H h NDy {0 o ) - z)- 1 + RT h (z) . (7.12) 

The right-hand side is the sum of three holomorphic terms in the interior of F h (t) and of an 
exponentially small term according to (|6.24l) . We obtain 

Po(t) = f (z- K)- 1 dz + O(e-t) = n (t) + 0(e"*) , 
znr J r h(i) 

where n (i) is the orthogonal spectral projector associated with {A^(t), . . . , A£(t)} C cr(if^(t)) 
with norm ||n (£)|| < 1 . 
d) For dtPo(t), we use 

- ffflo.v-wrH^o))- 1 = (z - H eo y_ wh (e ))-\d t V - d t W h ){z - H e0tV _ W H{e Q ))- 1 . 

From ([7TT2"]) . we get 

= + r h (z,V)l laM +RT h (z)l [aM 

4 3 

= E E y ) - A )il B ^ (l^z, 0~o), -) L2{R) l(e z ,z, 9 ) + RT h (z)l [aM , 

k=2 

where the first term of the right-hand side is holomophic inside T h {t) and the last term is exponen- 
tially small according to (|6.24[) and (|6 . 25|) . A symmetric writing holds for lr a ,6] — He .v-w h (^o))^ 1 ■ 
Hence the derivative dtPo(t) is the sum of several terms: 



1 



(z - K)-'(d t v - d t w n )(z - h*)- l dz = d t n (t) 

r"(t) 

' / {z-H^idtV -dtW^iz-H^)- 1 dz, (7.13) 

| Z -A°(t)|=c/2 



2l7T 

1 

2ri 



T h (z,V)l [a . b] (d t V - d t W h )(z - H h D ) dz 

r"(t) 



= -J2 ? h (X$,V)(d t V- d t W h )\^,){^,\ , (7.14) 
/ (z- J ff eo , y _^(0o))" 1 (a 4 y-9 t Ty' l )i?T' i (z) dz, (7.15) 



2i?r 



r*(t) 



plus another term symmetric to (|7.14l) . 

The first one (|7.13[) is uniformly bounded because 



• 11(2 ' 1 1 



Hd) 1 \\c(l 2 ) i s uniformly bounded when \z — A°(t)| = | according to Hypothe- 
sis canes) with HftVIU- + ||fitW?IU- < i , 

• ||(z — -Hd) -1 ^^ ^f l h ) is uniformly bounded when |z — A°(t)| = § according to (|7.4j)(|7.5D 
and (fA~T|) in Lemma [Q with \\dtW%\\ Mb < \ ■ 

The last one (I7TT5|) is 0{eT^k) owing to (|B^5]I for R$ h (z) and owing to ([g^T ]) for 

( z _ ■#e o ,v-iv fc (0o)) _ • 

For the middle term (|7.14p and its symmetric counterpart, first consider for k = 2, 3 and j' G 
{!,...£} 

{lie*, z, Oo) , (d t V - d t W h )^,) = ±^(u k , (d t V - d t W h )<$><],) , 
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where u-i and U3 are recalled in Lemma 16.41 The exponential decay estimates for stated in 
([6~T3))(|634)) and the one for stated in Proposition |4J}ii) combined with \\d t (V - W h )\\ Mb < \ 
imply that the scalar product is smaller than e~^h . Since the other factors of (|7.14[) are bounded 
by Ch or -M^, we conclude (I7.14[) and its symmetric counterpart are smaller than e~*£ . □ 



A Parameter dependent elliptic estimates on the interval 

[a, 6] 

We gather here elementary /i-dependent estimates for the elliptic operator — h 2 A + V on the interval 
(a, b). 

A.l Dirichlet problem 

It is convenient to use the /i-dependent 77 fc -norms 



2 _ 
fe , h — 

a<k 

for k € N. The estimates with the standard H k , k £ N, can be recovered after 



1 

For k = — 1, the ^.-dependent norm on i? _1 ((a, 6)) = [ZZq ((a, 6))] ' is 

■I, II !(/,")! 

\\j\\H- 1 . h ((a,b)) - S UP Tj— J] , 

u£Hi((a,b)) \\ U \\H^((a,b)) 

with now ||/|| ff -i,h(( a ,6)) < ^||/||iy-i(( tt ,6)) ■ We will note H^' h ((a,b)) the space H$((a,b)) equipped 
with the H 1 ' norm. 

Lemma A.l. There exists a constant C > such that 

C C 
Vw eC°°([0,h]), \u{0)\ < j-^\\u\\ H i, h{(0 ^ h)) , resp. \u' (0)\ < -^\\u\\ H 2,h((o,h)) > 

and the inequality extends to 1 ((0, /i)) (resp. H 2 ((0, h))) . 

Proof: The second estimate is simply a consequence of the first one after replacing u with hu' . 
The first estimate is simply the usual estimate \v(0)\ < C|| f || jyi ((0,1)) applied with v(x) = u{hx) . 

□ 

Lemma A. 2. Let Vi € L°°((a,b)) and V2 G At&((a, b)) be real valued with suppV^ CC (a, b) and 

IIVilUoo < i , \\v 2 \\ Mb <-. 

c c 

Then the Dirichlet Hamiltonian Hp = —h 2 A + Vi + V2 defined with the form domain Hq((cl, b)), 
satisfies the resolvent estimate 

Vz $ a{H h D ) , \\{H h D z)-X( H -^ ( (a m-Hb-Ua < Ca, b ,c [1 + \*\? U 



(A.l) 

When f € L 2 ((a, b)) and z £ er(i/^>) i/ie traces u'(a) and u'(b) of u = (Hp — z) 1 / are weZZ defined 
with 

\u'(a)\ + l«'0»l < (l + SjT^y) ll/lk. ■ (A.2) 
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Proof: From the Gagliardo-Nirenberg estimate |u(a;)| 2 < ^p-||nw'||£2||u|| L 2 with 1 1 V2 1 1 jVi fe 
the term (u , V2U) in the variational formulation is bounded by 

\{u, v 2 u)\<^\\hu'\\ L 4u\\ L 2<hhu'" 2 (r,J '" 



< 2. 



_ li^-r-^-w 
With 

Vu e Ht((a,b)) , (u,H*u) + C\\u\\h > \\\hu'\\%> + \ C - ||Vi||l- - N| 2 L 2 , 

the operator + C is bounded from below by -^-Aij + % when C > ~ + -§jr > ||Vi||i« + -§£ . 
Lax-Milgram theorem then says 

\\(H D + C) _1 /|| H -i.'»(( a , )i) )) < C a ,h iC ||/|| fl --i,h(( Q ,b)) • 

From the iterated first resolvent formula, 

(h h d - z)- 1 = (Hp + cy 1 + (c + z)(Hk + cy 2 + (c + zf{H h D + cy\H% - zy\H h D + cy 1 , 

we deduce (jA.ip . 

It contains also the estimate \\(Hp + Cy 1 \\ c( ^ L2 ^ a &)).#- 1 .'>(( a t))) — ~ Z h b " an< ^ w ith 
[H h D - zy 1 = (Hp + Cy 1 + (C + z){H h D + cy\H h D - 

this yields 

rh _x-i„ ^C aAc [l + \z\] 



\(H D z) 1 \\r(L 2 ((a.b)):Hb h ((a.b)Y) — 



When u — (Hp — z)^ 1 f with / G L 2 ((a, &)), writing the equation in [a, a + h] and [b — h, b] in the 
form —h 2 u" = / — (Vi — z)u implies 

11 h , C a .b c [l + \z\} 2 ( 1 \ 

lhlU=>((a,a +ft )U(6- M )) < ft ^ + d(z>ff (fl*)) J ll/lli2 ■ 

Lemma |A~T1 is applied to u(a + .) and u(b — .) in order to get (|A.2p . □ 
A. 2 Agmon estimate 

The next estimate is the usual energy estimate with exponential weights (see [2] [35]) . 

Lemma A.3. Let (a,/3) 6e an open interval, V € L°°((a,/3)) ; z e C and 95 e W 1,0 °((a, 0); K) . 
Denote by P the Schrddinger operator P := —h 2 d 2 /dx 2 + V. T/ien /or any u\,U2 in ff 1 ((a,/3)) 
sitc/i </ia< it" is a bounded measure in (a, /3) and locally L 2 around a and /3, the identity 

P_ rP 

U2e 2 h (P — z)u\dx = I hv^hv'idx + / (V — z — ip' 2 )v 2 vidx 

: J a J a 

y 

+ / hip' (v2v[ — v' 2 vi)dx 



+h 2 (e^^u'M) ~ e 2 ^U2uyp)) (A.3) 

holds by setting Vj := e ip / h Uj for j = 1, 2 . 

This identity is obtained after conjugation of hd/dx by e v ^ h and integration by parts. The 
weak regularity assumptions can be checked after regularizing individually u-y, U2, or V. In 46! 
it was even considered with possible jumps of the derivative at a and /3, which are here removed 
by the simplifying condition that u\ is locally H 2 around a and (3 (Jump conditions already occur 
at the ends of our intervals). 
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B Variation on adiabatic evolutions. 

We shall consider a family of contraction semigroup generators (A e (s)) s6 [ 0i +oo) which fulfill the 
two next properties. 

• The Cauchy problem 

ied t u t = iA £ (t)u t ,^ ^ 

u t= o = mo 

admits a unique strong solution with u t G D(A £ (t)) for all t > as soon as wo € -D(^ e (0))- 
The corresponding dynamical system of contractions is denoted (S £ (t,s)) t > s with the prop- 
erty S £ (t,s)D(A £ {s)) c D{A £ {t)). 

• The resolvent (z - iA e (s))~ 1 defines C x+1 ([0, +00); £(%)) function for some z G C and that 
the exists a contour T £ C C independent of s G [0,T], such that 

|r e |+ max jia^-^))- 1 !!^^, 

for any k G {0, . . . , K + 1}, K G N, any (5 G (0, <5 ) and any e G (0, e ) . 
Notation: We shall use the notation 17(e) = 0(e N ) for any N G Z in order to summarize 

V<5 >0,3C S;5 >0, | ff (e)| <C B , s e N - 5 . 
For example, the previous assumption can be written 

|r e |=0( £ °) and max Jd^z - ^(s^W = O(e ) . (B.2) 

The spectral projection Po(t) = Eo(t) is defined as a contour integral along T £ of the resolvent 
[z — A £ (t))~ 1 . Correction terms Ej(t), 1 < j < K are then constructed by induction. The finite 
sequence (E £ ) <j<K is defined according to 

E s (s) = P§{*) = ^J r ( z - ^W)" 1 rfz > OoW = 1 - Po(s) , (B.3) 
S £ (s)=J2 E£ m(s)E £ j - m (s), i£2< j<K, SZ = St = 0, (B.4) 

m— 1 

E s j(s) = ^J ^ {QoW^-i(*)^W - P £ («)^^-i(«)Qo(«)} ^ <** (B.5) 

- 2P 5 (s)5|(s)P £ (s) , with i? e = (z - A e (*)) _1 • (B.6) 

Theorem B.l. There exists a C -projection valued function (P £ (s)) s >q such that the relations 
and estimates 

P £ (s)P £ (s) = P £ (s) , P £ (s)eC(n,D(A £ (s))), (B.7) 
if 

P £ {s) = ^e J E £ (s)+d(e K+1 ) with E £ (s) = O(e ) m £(«) , (B.8) 

P £ (s)A £ (s) =O(e ) and ^ £ (s)P £ (s) = O(e ) in £.(U) , (B.9) 
ied s P £ (s) - [iA e (s), P £ (s)} = d(e K+1 ) in C(H) , (B.10) 

ftoZd with uniform constants with respect to s,t G [0,T] /or any /ixed T < +00 . Moreover for 
K > 1, if u s = P £ (s)u s then u £ (t) = S £ {t 1 s)u s satisfies 

sup \\u £ (t)-v £ (t)\\=0(s K ), (B.ll) 

s<t<T 

with v £ (t) = P £ (t)v £ (t), (B.12) 
and ( iftV s -ie(d t P £ (t))v £ = P £ (t)(iA £ (t))P £ (t)v £ , for t > s , (Rig) 
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The proof of this theorem follows the lines of [33]. For the sake of completeness, we check that the 
computations are still valid in the non self-adjoint unbounded case (bounded self-adjoint generators 
have been considered in [44] [55]) and that the O(e ) estimates can be propagated in the induction 
process like the uniform constants in [33] . Part of the analysis could be pushed further in the spirit 
of [35] in order to get (D(e~i) error under analyticity assumptions but the techniques developed by 
A. Joye in this article should be adapted in order to work with a e-dependent gap or with O(e ) 
resolvent estimates, maybe by including all the additional information provided by our model. 

As it is stated, the previous result cannot be used for K — and is not formulated as usual 
with a reduced evolution on the fixed space RaixP e (s) after introducing the parallel transport 
associated with the C 1 family (P £ (t))t> s ■ Actually both problems can be solved at the first order 
with an additional uniform boundedness assumption on Eq{£) and dtE^t) . This will be obtained 
as a corollary of Theorem lB.il used with K = 1 before reconsidering the case K = 0. The parallel 
transport s'), associated with (-Eo(i))te[o,Tl) is defined for t',s' G [0,T] by 

\ #§(f = *V)=M, ( ' 

and the uniform boundedness of $5(t' , s') is inherited from the one of E^t) and &tE^{t) . 

Corollary B.2. With the hypotheses of Theorem \B.1\ with K > 1, assume additionally that the 
projector Eq(s) defined in (|B.3I) and its derivative <9 s _Eq(s) are uniformly bounded continuous func- 
tions: 

3C> 0,Ve e (0,e ), max \\E^s)\\ + \\d s E e (s)\\ < C . 

s£ [0,T] 

Then for K > 1 and when u s = E^(s)u s , the solution u e {t) = S e (t, s)u s to (|B.1[) satisfies 

sup \\u 6 (t) - $g(t, s)w e (t)\\ = d(e)\\u s \\ , (B.15) 

s<t<T 

where $o(^'j s ') * s ^ e parallel transport defined for t 1 , s' G [0,T] by (|B.14[) and w £ € Eq(s)'H solves 
the Cauchy problem 

ied t w £ = %(s,t)E^t)(iAHt))E £ (t)%(t,s)w^ = £J§(«)*8(M)(^ e (*))*§(*, s)E s (s)w s 
w e (t = s) = u s . 



Theorem IB . 1 1 and Corollary IB.2I are proved in several steps. We start with uniform estimates 
for the Ej's. 

Proposition B.3. For all j 6 {0, . . . , K}, and any T E M + , the £(W)-valued functions Ej and 
Sj satisfy: 

K+l-j 

\\d*EZ(s)\\ + \\d*S^ S )\\=O(e ), E 3 (s)eC(n,D(A^s))), (B.16) 

fc=0 

A E (s)E j (s) = O(e ) and E J (s)A E (s) = 5(e°) in C(H) , (B.17) 
E£ M = E Et n (s)EZ_ m ( s ) = S*(s) + EUsW^s) + E$(s)E%(a) , (B.18) 

m— 

id.Ej^is) = [iA s (s), , for j > 1 , (B.19) 

with uniform constants w.r.t. s € [0, T] . 

Proof: The first statement for j — is a consequence of the definition (|B.3j) of Eq(s) = Pq(s) 
combined with the estimates (|B.2I) of d^(z — iA £ (s))~ 1 . By induction assume that the properties 
are satisfied for j < J < K . The definition (|B.4|) of Sj+i an< ^ (jB.5[) of £Jj , 1 provide directly the 
first statement (|B.16[) for j = J + 1 . The second statement of (|B.16[) and the estimates (|B.17I) 
rely on the bound of the J r -term which is obtained after noticing 

A^SW = -1 + {z _ Z A e {8)) = RCA ^ S ) ( B ' 2 °) 
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and the bound of the two other terms which is deduced from the induction assumption for j < J . 
Compute the commutator [iA £ (s), Sj +1 (s)]: 

J 

[iA s , S e J+1 ] = J2 (K> E e J+1 _ m ] + [iA s , E^} P} +1 _ m ) , 

m— 1 
J 

= Y, E Uid s E E j_J + {id a E^ l )E l j +l _ m , owing to fEU} with j < J , 

= id s E E j - E e Q {id s E^j) - {id s E e j)E e , owing to (|B~T8j) with j < J , 
= -i (Pi(d s Efj)P^ - Ql(d a E e j)Ql) , owing to E s = P e = 1 - Q% . 

With P e P e = P £ and P e Q§ = QgP e = 0, this implies 

[iA s , P e 55 +1 P £ ] = ~iP*(d s E^)Pl , (B.21) 
[iA s , Q e S £ J+ iQo] = iQl{d s E s j)Ql , (B.22) 
[iA°, PiS s J+1 Ql] = [iA s , QIS} +1 P%] = . (B.23) 

The definition of Pq as a spectral projection associated with iA £ and (|B.23|) . imply (see for example 
[33] Proposition 1) 

PoSj+iQq = QqSj+xPq = • 
Meanwhile the definition (|B.5|) of Ef J+1 for j = J + 1 implies 

^0 ^,7+1-^0 = — ^O^J+l-^O > QqEj + iQo = QqSj +1 Qq . 

This yields the relation (|B.18|) for j = J + 1. Another consequence with (|B.21|) and (|B.22|) is 

[iA°, P £ P5+i^o + QoE e J+ iQo] = P^idsE^PI + Ql(idE £ j)Ql . (B.24) 

Finally compute the off-diagonal blocks Pq [iA e , Ej +1 ]Qq and Qq[iA s , E £ j+1 ]Pq by using again the 
definition ([53]) of E e J+1 , the relation (fIQ3"1) . and the identity ((H7201) : 

[iA s , i^+iQo] = P E [iA s , E S J+1 ] Q% = »ifi(0.£5)Q§ , 
[tA £ ,Q £ Q E £ J+1 P £ ] = iQKd.E°j)P§ . 

Summing these last two equalities with (IB.24[) yields the relation (|B.19I) for j = J + 1 . □ 
The above calculations are essentially the same as in [43 33j and, as a consequence of Proposition 
IB. 31 the sum 

K 

T s («) = ^Vl^) (B.25) 

solves (|B.7|) . (|B.9|) and (|B. 10[) in the sense of asymptotic expensions; in particular: 

T e 2 = T e + d(e K+1 ) (B.26) 

Here comes the main difference which is necessary because no better estimate than ||Pq || = O(e ) 
can be expected in our non self-adjoint case. We will need the next lemma. 

Lemma B.4. Assume that T € C(H) satisfies \\T 2 - T\\ < 5 < 1/4 and \\T\\ < C then 

a{T) C {z e C , \z(z - 1)| < S} C {z £ C , \z\ < c s } U {z e C , \z - 1| < c s } , c s < 
max|||(z-T)- 1 ||,|z-l| = i|<2^ 



1 

2 ' 
2 C+1 

^45' 



2(7 +1 If 
|r-P||<2 - , f 3, iwtfi P=— / fjs-T) _1 dz 
1-45 2t7r ./|,-i|=i/2 
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Proof: If z G a{T) then z(z - 1) belongs to a(T(T - 1)) c{z£C, \z(z - 1)\ < \} (Remember 
that |z(z — 1)| = t means \Z — 4| = \ with Z = (z — \) 2 ). Consider z & C such that \z — 1| = \, 
then the relation 

(T-z)(T-(l-z)) = z(l-z) + (T 2 -T), 
with \z(l -z)\>\ and \\T 2 - T|| < S < \, implies 

II (T - z)- 1 !! < ||T - (1 - z)||||[z(l - z) + T 2 - T]- 1 !! < for \z 1| = 1 . 

The symmetry with respect to z = | due to (1 - T)(l - T) - (1 - T) = T 2 - T implies also 

\\{T-z)- 1 \\< ( ^\ for |*| = £. 



4 * 



Compute 

P = 

2in 



T-P = ±-f [T(z-l)-i-(z-T)-i] dz 

|a-l|=* 



= {T -1)P +(T 2 -T)A X with A 1 = —[ (T - z)~ l (l- z)~ l dz . 

In particular this implies to T(l - P) = (T 2 - T)Ai while replacing T with (1 - T) and P with 
1 — P leads to 

P - T = -T(l - P) + (T 2 - T)A with A = — [ {T - z)- 1 z- 1 dz . 

2t7r J\ z \=i 

We finally obtain 

T -P =(T 2 - T)(A 1 - A ) = {A! - A )(T 2 - T) (B.27) 
and ||T - P|| < \\T 2 - T\\ [||Ai|| + ||A ||] < 5 x . 



4 

□ 

Proposition B.5. Consider the approximate projection T s (s) defined in (|B.25I) . then there exists 
a projection P £ {s) such that 

\\P £ (s) - T e (s)\\ = d(e K+1 ) , P £ (s) g £(H, D(A%s))) , (B.28) 
\\A e (s)P e (s)\\=d(£°) and \\P £ (s)A £ (s)\\ = O(e ) , (B.29) 
ied s P £ (s) = [iA £ (s), P £ (s)} + 0{e K+1 ) , m C(H) . (B.30) 

Proof: For e > small enough, set 

P e (s) = ^- / (z-T^))" 1 ^. 

Owing to (IB.26|) . the first statement of (|B.28|) is a straightforward application of Lemma TB. 41 The 
definition of T s (s) implies A £ (s)T e (s) G C(U) and T £ (s)A £ (s) G £(%). The relation (|B~27j) with 
T = T £ and P = P £ gives (f5T2^|> . Computing the derivative d s P £ {s) with T e G C 1 ([0, T]; £(%)) 
gives: 



ie^P e = (z-T £ )- 1 (ze5 s r e )(z-r e )- 1 dz 

= ^ I (z ~ Te)- 1 [iA £ , T e ] (z - T,)- 1 dz 

+ / — ^Ir) (ds-EfcX* — T e ) dz . 
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The relation 

[iA £ , (z - T,)- 1 ] =(z- T,)- 1 [iA £ ,T e ] (z ~ T,)- 1 

allows to conclude 

ied s P £ - [iA*,P°] = 4t— / (z - T £ )~ x (d s E e K )(z - T,)- 1 dz = 0(e K+1 ) . 

□ 

Proof of Theorem IB.lt The statements (jB.7p(jB.8p and (jB.lOp have already been checked in 
Proposition IB. 31 and Proposition IB. 51 Consider now the adiabatic evolution of S(t,s)u s , when 
u s = P e (s)u s , stated in (IBTTj) (jB~T2j) (jB~T3)) . 

We assumed that the Cauchy problem (|B.1[) defines a strongly continuous dynamical system 
(S e (t,s))t> s >o of contractions in % with S e (t, s)D(A £ (s)) C D(A £ {t)). We now consider the 
modified operator 

H% D (t)=iA%t)+B%t), 
with B £ (t) = (1 - 2P £ (t)) (ied t P £ (t) - [iA e (t), P £ (t)}) . 

Since B £ (s) is an 0(e K+1 ) bounded continuous perturbation of iA £ (s), the Cauchy problem 

f ied t ut = (iA £ (t)-B £ (t))u t 

\ U t = s = U s , 

defines a strongly continuous dynamical system of bounded operators 

S AD (t, s) = S £ (t, s) - is" 1 J S £ (t, s>)(-B%s'))S AD (s', s) ds' . 

With the help of Gronwall lemma, it satisfies 

\\S e AD {t,s)\\ < e' -1 /.' H^COII ^' < e C s , T e«S{t-s) 
\\S £ (t, s) - S £ AD (t, s)\\ < C StT e K -\t - s)e c ^ £K "^-^ , 
for all s, t, < s < t < T, with S £ (0, 1). Note that the right-hand side is bounded when 

K > 1 . 

For the comparison (jB.111) we take simply v £ {t) = S £ AD (t, s)u s . It remains to check (|B.12j) and 
(IB.13P - First notice the identity 

H% D (t) = P e (t)(iA s (t))P e (t) + (l-P%t))(iA^t))(l-P e (t)) (B.31) 
-ie [P £ (t)d t P £ (t) + (1 - P £ (t))d t (l - P £ (t))} . 

For our choice v £ (t) = S £ AD (t,s)u s the quantity P £ (t){ied t v £ (t)) equals P e (i)(iA e (i) + B £ {t))v £ (t) 
since P e (i)A e = A £ (i)P £ (i) - [A £ (t), P £ {t)] e £(%) . Hence P £ (t)v £ (t) satisfies in the strong sense 
the equation 

(ied t -H £ AD (t))(P £ v £ ) = ie(d t P £ )v £ +P £ (H £ AD (t)v £ )~H £ AD (t)P £ (t)v £ 

= ie [{d t P £ ) - P £ {d t P £ ) - {d t P £ )P £ ] v £ {t) = . 

As a strong solution to idtv = H AD (t)v(t) with the initial data P £ (s)u s = u s , P £ (t)v £ (t) has to be 
equal to v £ (t) . We have proved (IB. 121) . The equation (|B.13|) is a rewriting of iedtv £ — H A o (t)v £ = 
after recalling n(<9II)n = when If 2 = II . □ 
Proof of Corollary IB.2t Theorem IB. II applied with K = 1 gives the approximation v £ (t) = 
P £ (t)v £ (t), P £ {t) = E £ (t) + eEf(t) + die 2 ), which solves 

f ied t v £ - ie(d t P £ (t))v £ = P £ (t)(iA £ {t))P £ (t)v £ , for t > s , 
\ v £ (t = s) = u s . 
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The relations P £ = E £ + eE{ + C(e 2 ), [A £ , E £ ] = and E £ E(E £ = combined with the estimates 
(|B~T71) and (|R29|l lead to: 

P £ mA £ (t))P £ (t) = E £ mA £ (t))E £ (t) + eEimA £ {t))El{t) + e{iA £ {t))E £ {t)El{t) + d(e 2 ) . 

Then Proposition IB. 31 provide 

ied t P £ (t) = ied t E £ (t) + d(e 2 ) . 

This implies 

ied t v £ ~ ie(d t E £ {t))v £ = E £ (t)(iA £ (t))E £ {t)v £ + eEf (t)(iA £ {t)E £ {t))v £ + d(e 2 ) , (B.32) 
where we used E £ v £ = v £ + 0(e). Consider now the adiabatic generator 

H £ {t)=iA £ {t)+B £ {t), 
with B £ (t) = (1 - 2E £ (t)) (ied t E £ (t) - [iA £ (i), E £ (t)\) = ie(l - 2E £ (t))d t E £ (t) . 

The assumed estimates on Eq and dtEg with the Gronwall Lemma lead to the uniform bound for 
the associated dynamical system S^it, s): 

Vs,t,0<s<t<T, \\S £ (t,s)\\ <C T e CT , 

while the formula (|B.31| is valid for H £ {t) after replacing P £ (t) with E £ (t): 

H £ (t) = E £ mA £ (t))E £ (t) + (1 - E*(t))(iA £ (t))(l - E £ (t)) 

- l£ [E £ (t)d t E £ (t) + (1 - E £ (t))d t (l - E £ (t))} . 

Now compute 

(ied t - H £ (t))(E £ v £ ) = ie{d t E £ )v £ + E £ (ied t v £ ) - H £ (t)E (t)v £ . 

With (|ET32|t and E^E £ E £ = 0, one gets 

{iedt - H £ (t))(E £ v £ ) = -is [(d t E £ ) E £ {d t E £ ) ~ (d t E £ )E £ ] v £ (t) + d(e 2 ) = d(e 2 ) . 

The uniform estimate ||<So(f, s )ll — Cye implies 

\\u £ (t)-S £ (t,s)u s \\ < \\u £ (t)~v £ (t)\\ + \\eE £ 1 (t)v £ (t)\\ + \\E £ (t)v £ - S £ (t, s)u s \\ = 5(e) . 

The same argument as in the end of the proof of Theorem IB. II says that = S £ (t,s)u s with 

E £ (s)u s — u s satisfies 

Vie[ S) T], E £ (t)v £ (t) = v £ (t) 

and solves the Cauchy problem 

f ied t v £ - ie(d t E £ (t))v £ Q (t) = E £ (t)(iA £ (t))E (t)v £ (t) , 
\ v^(t = s) = u s . 

The uniform boundedness of E £ (t) and dtE^t) ensures that the solution to (IB. 14[) is well defined 
for t', s' € [0, T] with the uniform estimate 

vt' )S 'e[o,T], \\m',s')\\<C' T , 

with the parallel transport property 

W,s'e[0,T], mt', S ')E s ( S ')=E*(t')mt',s'), [^y)]" 1 = #§(«',«'). 

It suffices to take w £ (t) = <&q(s, t)S £ (t, s)u s . □ 
Acknowledgements: The authors were partly supported by the ANR-project Quatrain led 
by F. Mehats. This work was initiated while the second author had a CNRS post-doc posi- 
tion in Rennes. The third author discussed recently this problem with G.M. Graf, E.B. Davies, 
K. Pankrashkin and he remembers many other former discussions related to this topic with col- 
leagues appearing in the bibliography. 



56 



References 



W.K. Abou Salem and J. Frohlich. Adiabatic theorems for quantum resonances. Comm. 
Math. Phys., 273(3):651-675, 2007. 

S. Agmon. Lectures on exponential decay of solutions of second-order elliptic equations: bounds 
on eigenf unctions of N -body Schrodinger operators, volume 29 of Mathematical Notes. Prince- 
ton University Press, Princeton, NJ, 1982. 

J. Aguilar and J. M. Combes. A class of analytic perturbations for one-body Schrodinger 
Hamiltonians. Comm. Math. Phys., 22:269-279, 1971. 

S. Albeverio and K. Pankrashkin. A remark on Krein's resolvent formula and boundary 
conditions. J. Phys. A, 38(22) :4859-4864, 2005. 

X. Antoine, A. Arnold, C. Besse, M. Ehrhardt, and A. Schadle. A review of transparent 
and artificial boundary conditions techniques for linear and nonlinear Schrodinger equations. 
Commun. Comput. Phys., 4(4):729-796, 2008. 

A. Arnold. Mathematical concepts of open quantum boundary conditions. Transp. Theory 
Stat. Phys., 30(4-6):561-584, 2001. 

J. E. Avron, A. Elgart, G. M. Graf, and L. Sadun. Transport and dissipation in quantum 
pumps. J. Statist. Phys., 1 16(1-4) :425-473, 2004. 

J. E. Avron, R. Seiler, and L. G. Yaffe. Adiabatic theorems and applications to the quantum 
Hall effect. Comm. Math. Phys., 110(l):33-49, 1987. 

J.E. Avron, A. Elgart, G.M. Graf, L. Sadun, and K. Schnee. Adiabatic charge pumping in 
open quantum systems. Comm. Pure Appl. Math., 57(4):528-561, 2004. 

E. Balslev and J. M. Combes. Spectral properties of many-body Schrodinger operators with 
dilatation-analytic interactions. Comm. Math. Phys., 22:280-294, 1971. 

M. Baro, H.C. Kaiser, H. Neidhardt, and J. Rehberg. Dissipative Schrodinger-Poisson systems. 
J. Math. Phys., 45(l):21-43, 2004. 

M. Baro, H.C. Kaiser, H. Neidhardt, and J. Rehberg. A quantum transmitting Schrodinger- 
Poisson system. Rev. Math. Phys., 16(3):281-330, 2004. 

M. Baro, H. Neidhardt, and J. Rehberg. Current coupling of drift-diffusion models and 
Schrodinger-Poisson systems: dissipative hybrid models. SI AM J. Math. Anal, 37(3) :941- 
981, 2005. 

J. Behrndt, M.M. Malamud, and H. Neidhardt. Scattering theory for open quantum systems 
with finite rank coupling. Math. Phys. Anal. Geom., 10(4):313-358, 2007. 

N. Ben Abdallah. On a multidimensional Schrodinger-Poisson scattering model for semicon- 
ductors. J. Math. Phys., 41(7):4241-4261, 2000. 

N. Ben Abdallah, P. Degond, and P.A. Markowich. On a one-dimensional Schrodinger-Poisson 
scattering model. Z. Angew. Math. Phys., 48(1):135-155, 1997. 

N. Ben Abdallah and O. Pinaud. Multiscalc simulation of transport in an open quantum 
system: resonances and WKB interpolation. J. Comput. Phys., 213(1):288-310, 2006. 

V. Bonnaillie-Noel, A. Faraj, and F. Nier. Simulation of resonant tunneling heterostructures: 
numerical comparison of a complete Schrodinger-Poisson system and a reduced nonlinear 
model. J. Comput. Elec, 8(1):11-18, 2009. 



57 



[19] V. Bonnaillie-Noel, F. Nier, and Y. Patel. Computing the steady states for an asymptotic 
model of quantum transport in resonant heterostructures. J. Comput. Phys., 219(2):644-670, 
2006. 

[20] V. Bonnaillie-Noel, F. Nier, and Y. Patel. Far from equilibrium steady states of 1D- 
Schrodinger-Poisson systems with quantum wells. I. Ann. Inst. H. Poincare Anal. Non 
Lineaire, 25(5):937-968, 2008. 

[21] V. Bonnaillie-Noel, F. Nier, and Y. Patel. Far from equilibrium steady states of 1D- 
Schrodinger-Poisson systems with quantum wells. II. J. Math. Soc. Japan, 61(1):65-106, 
2009. 

[22] H. D. Corncan, P. Duclos, G. Nenciu, and R. Purice. Adiabatically switched-on electrical bias 
and the Landauer-Biittiker formula. J. Math. Phys., 49(10):102106, 20, 2008. 

[23] H. L. Cycon, R. G. Froese, W. Kirsch, and B. Simon. Schrddinger operators with application to 
quantum mechanics and global geometry. Texts and Monographs in Physics. Springer- Verlag, 
Berlin, study edition, 1987. 

[24] M. Dimassi and J. Sjostrand. Spectral asymptotics in the semi- classical limit, volume 268 of 
London Mathematical Society Lecture Note Series. Cambridge University Press, Cambridge, 
1999. 

[25] M. Ehrhardt and A. Arnold. Discrete transparent boundary conditions for the Schrodinger 
equation. Riv. Mat. Univ. Parma (6), 4*:57-108, 2001. Fluid dynamic processes with inelastic 
interactions at the molecular scale (Torino, 2000). 

[26] H.O. Fattorini. The Cauchy problem, volume 18 of Encyclopedia of Mathematics and its 
Applications. Addison- Wesley Publishing Co., Reading, Mass., 1983. 

[27] C. Gerard and I. M. Sigal. Space-time picture of semiclassical resonances. Comm. Math. 
Phys., 145(2):281-328, 1992. 

[28] B. Hclffcr. Semi- classical analysis for the Schrodinger operator and applications, volume 1336 
of Lecture Notes in Mathematics. Springer- Verlag, Berlin, 1988. 

[29] B. Hclffcr and A. Martinez. Comparaison entre les diverses notions de resonances. Helv. Phys. 
Acta, 60(8):992-1003, 1987. 

[30] B. Hclffcr and J. Sjostrand. Resonances en limite semi-classique. Number 24-25 in Mem. Soc. 
Math. France (N.S.). 1986. 

[31] P. D. Hislop and I. M. Sigal. Semiclassical theory of shape resonances in quantum mechanics, 
volume 78(399) of Mem. Amer. Math. Soc. 1989. 

[32] P. D. Hislop and I. M. Sigal. Introduction to spectral theory, volume 113 of Applied Mathemat- 
ical Sciences. Springer- Verlag, New York, 1996. With applications to Schrodinger operators. 

[33] G. Jona-Lasinio, C. Presilla, and J. Sjostrand. On Schrodinger equations with concentrated 
nonlinearities. Ann. Physics, 240(1):1-21, 1995. 

[34] A. Joye and C.E. Pfister. Exponential estimates in adiabatic quantum evolution. In Xllth 
International Congress of Mathematical Physics (ICMP '97) (Brisbane), pages 309-315. Int. 
Press, Cambridge, MA, 1999. 

[35] Alain Joye. General adiabatic evolution with a gap condition. Comm. Math. Phys., 275(1):139- 
162, 2007. 

[36] T. Kato. Integration of the equation of evolution in a Banach space. J. Math. Soc. Japan, 
5:208-234, 1953. 



58 



[37] T. Kato. Linear evolution equations of "hyperbolic" type. J. Fac. Sci. Univ. Tokyo Sect. I, 
17:241-258, 1970. 

[38] M. Klein and Rama J. Almost exponential decay of quantum resonance states and Paley- 
Wicner type estimates in Gevrey spaces, preprint, mp-arc 09-64, 2009. 

[39] M. Klein, J. Rama, and R. Wiist. Time evolution of quantum resonance states. Asymptot. 
Anal, 51(1):1-16, 2007. 

[40] A. Lahmar-Benbernou and A. Martinez. On Hclffer-Sjostrand's theory of resonances. Int. 
Math. Res. Not, (13):697-717, 2002. 

[41] S.E. Laux, A. Kumar, and M.V. Fischetti. Analysis of quantum ballistic electron transport in 
ultra-small semiconductor devices including space-charge effects. J. Appl. Phys., (95):5545- 
5582, 2004. 

[42] P. Lochak. About the adiabatic stability of resonant states. Ann. Inst. H. Poincare Sect. A 
(N.S.), 39(2):119-143, 1983. 

[43] G. Nenciu. Linear adiabatic theory. Exponential estimates. Comm. Math. Phys., 152(3):479- 
496, 1993. 

[44] G. Nenciu and G. Rasche. On the adiabatic theorem for nonselfadjoint Hamiltonians. J. Phys. 
A, 25(21):5741-5751, 1992. 

[45] F. Nier. The dynamics of some quantum open systems with short-range nonlinearitics. Non- 
Uneanty, 11(4):1127 1172, 1998. 

[46] F. Nier. Accurate WKB approximation for a ID problem with low regularity. Serdica Math. 
J., 34(1):113-126, 2008. 

[47] K. Pankrashkin. Resolvents of self-adjoint extensions with mixed boundary conditions. Rep. 
Math. Phys., 58(2):207-221, 2006. 

[48] G. Perelman. Evolution of adiabatically perturbed resonant states. Asymptot. Anal, 22(3- 
4): 177-203, 2000. 

[49] Olivier Pinaud. Transient simulations of a resonant tunneling diode. J. App. Phys., 
92(4):1987-1994, 2002. 

[50] C. Presilla and J. Sjostrand. Transport properties in resonant tunneling heterostructures. J. 
Math. Phys., 37(10):4816-4844, 1996. 

[51] M. Reed and B. Simon. Methods of modern mathematical physics. II. Fourier analysis, self- 
adjointness. Academic Press, New York, 1975. 

[52] M. Reed and B. Simon. Methods of modern mathematical physics. IV. Analysis of operators. 
Academic Press, New York, 1978. 

[53] B. Simon. Resonances and complex scaling: a rigorous overview. Int. J. Quantum Chem., 
14(4):529-542, 1978. 

[54] B. Simon. The definition of molecular resonance curves by the method of exterior complex 
scaling. Phys. Lett., 71A(2,3):211-214, 1979. 

[55] J. Sjostrand. Projecteurs adiabatiques du point de vue pseudodiffcrentiel. C. R. Acad. Sci. 
Paris Ser. I Math., 317(2):217-220, 1993. 

[56] J. Sjostrand and M. Zworski. Complex scaling and the distribution of scattering poles. J. 
Amer. Math. Soc, 4(4):729-769, 1991. 



59 



[57] E. Skibsted. Truncated Gamow functions, a-decay and the exponential law. Comm. Math. 
Phys., 104(4):591-604, 1986. 

[58] E. Skibsted. Truncated Gamow functions and the exponential decay law. Ann. Inst. H. 
Poincare Phys. Theor., 46(2): 131-153, 1987. 

[59] E. Skibsted. On the evolution of resonance states. J. Math. Anal. AppL, 141(l):27-48, 1989. 

[60] A. Soffcr and M. I. Weinstein. Time dependent resonance theory. Geom. Fund. Anal., 
8(6):1086-1128, 1998. 

[61] K. Yosida. Functional analysis. Classics in Mathematics. Springer- Verlag, Berlin, 1995. 
Reprint of the sixth (1980) edition. 



60 



